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INTRODUCTION 


The cross section of a forest tree tells its life story. Annual rings 
of medium width near the center reveal the tree’s origin in an isolated 
position or an even-aged stand, for two-storied or many-aged stand 
conditions are shown by narrow ring widths. A series of wide rings 
usually follows, denoting the tree’s vigorous prime, and narrowing 
rings near the circumference indicate its senescence. Within these 
broad age periods, cycles of narrower or wider rings suggest epochs 
of suppression by neighboring trees or dominance over them in the 
intense struggle for light and moisture; or they may imply periodic 
fluctuations of climate. Within such epochs, individual annual rings 
commonly show striking variation in width around the mean width 
for the epoch. Typical fluctuations of annual ring width are pre- 
sented diagrammatically in figure 1. 
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FIGURE 1.—Graphical representation of fluctuations in width of annual rings of a tree. 


Insofar as competition with neighboring trees and other vegetation 
affects growth within an epoch, the magnitude of the growth of a 
given tree in an unthinned stand may be expected to be quite inde- 
pendent of that of other trees, chosen at random, during the same 
span of years; for two or more trees are probably not subject to the 
same hazards from this source during the same time intervals. All 
the trees of a stand, however, are exposed to the same general condi- 
tions of weather and periodic fluctuations of climate. Insofar as 
weather affects growth, therefore, the width of the annual rings of 
single trees, or the average width of those of groups of trees, should, 
in the same years, be highly correlated, particularly if the effects of 
the age cycle and of epochs of suppressive competition be eliminated. 

! Received for publication July 22, 1936; issued March 1937. 


? During most of a year’s leave of absence from Switzerland, Dr. Meyer acted as collaborator in the 
section of forest measurements of the Forest Service. 
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It has been generally recognized that variation in width of annual 
rings is associated with variation in climate (1, 3, 4, 5, 8, 9, 10, 12, 18). 
So far as the present authors know, however, no attempt has been 
made to investigate the magnitude of the variation in annual growth 
that is due to fluctuation in climate. Such a measure is basic to 
judgment concerning the efficacy of a given silvicultural treatment in 
a regulated forest; for the comparison of periodic growth before and 
after treatment may be misleading unless the forester has some knowl- 
edge of the extent to which fluctuation in weather and climate may 
influence growth without change of treatment. 

Meteorological conditions as a whole are among the most complex 
phenomena of nature, subject as they are to an infinite number of 
combinations, each of which may influence growth in a different way. 
While it would hardly be possible to explain their effects fully—that is 
to find a functional relation between weather and growth—part of 
the effect of weather on growth may perhaps be explained, not as a 
functional relation but as a probable or conjectural relation, by the 
method of least squares. It is desirable to know, for instance, in 
what seasons and to what extent tree growth is affected by rainfall, 
by temperature, by bright sunshine, and so on. Positive answers to 
such questions should form valuable background for work in timber- 
tree physiology and ecology and in the practice of silviculture. 

Meteorological phenomena cannot be dealt with in their entirety, 
but the effect of single factors on growth can be investigated. Usually 
any one such factor is correlated with others, which are, therefore, 
also dealt with, although indirectly; rainfall, for instance, is correlated 
with temperature and sunshine, and with relative humidity. 

In this paper it is proposed (1) to divide the variation of annual 
growth of sample trees over a period of 63 years into that portion due 
to influence within and that to influence beyond the stand, to arrive 
at the magnitude of annual growth fluctuations that are due to out- 
side influences, and to illustrate the resulting effects on changes in 
basal area and volume growth; and (2) to investigate in detail the 
effects of climate on timber growth, with special reference to rainfall 
distribution and its influence on annual-ring width. 


TREE DATA 


Dr. Meyer brought with him from Switzerland the tree data used, 
which were made available to him through Dr. Hermann Knuchel, 
professor of forestry at the Eidgenéssische Technische Hochschule, 
Zurich. They consist of measurements of annual rings of 12 domi- 
nant white firs (Abies alba Mill.) taken from an even-aged stand near 
Aarau, Switzerland, for the years 1864-1926. Although the species 
and conditions may not be entirely comparable to any American 
species or environment, the methods used are equally appropriate to 
American timber types. 

Concerning the trees and their environment, a translation of 
Knuchel’s statement (8, pp. 268-269) is as follows: 

a large number of dominant sample trees had to be cut in an approximately 120- 
year-old fir-spruce-beech stand of the Swiss central highlands. 

The stand in question forms part of the so-called ‘‘lower forest”’ of the Zofingen 


Commune, 480 meters above sea level on a gently inclined northern slope of a 
Tertiary sandstone hill. The marly-loam soil is principally of glacial origin. The 


§ Reference is made by number (italic) to Literature Cited, p. 107. 
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average annual precipitation based on a 50-year record is about 1,050 mm; the 
average annual temperature is 8.3° C. 

The stand had been only slightly thinned up to a few years prior to the removal 
of the sample stems; at intervals of 5 and 6 years suppressed and dying trees only 
had been removed. In the last 20 years the thinnings were somewhat heavier, 
but the stand in 1926 was still fairly dense. The rate of growth of the spruce had 
fallen off somewhat, whereas the firs were still in full development. This fact 
induced us for the present to use only Abies pectinata D.C. [Abies alba Mill.] 
in the investigation of increment fluctuations. 

Of the 12 sample stems cross sections were taken at ground level and at heights 
of 4,8, and 12m. The height of the trees was 32 to 38 m; the ciameter breast 
high was 54 to 68 em. The central rings showed that the firs had originated 
naturally and that cutting of the old stand had taken place when the firs were 
10 to 30 years old and 1 to 6m high. On the basis of the annual rings we are 
justified in assuming that the trees had subsequently always been dominant. 


Width of individual rings at each crosscut were measured along 
three radii removed from each other by about 120°, and the average 
of these for each year entered as the tree-year observation. The 
present data are the measurements taken on the crosscuts 4 m above 
the ground. These are presented graphically in figure 2. Even in 
this raw form, correlation between trees in the different years is 
apparent. Deviations of annual ring widths from trends of the 12 
individual trees for the 63-year period are given in table 1, the last 
column of which lists all the variation that is due to meteorological 
conditions from year to year plus an element due to random error. 

The rainfall data are the monthly precipitation records at Aarau, 
Switzerland, from June 1863 to August 1926 (14). 


TABLE 1.—Deviations of annual ring width from trends of individual trees, in 
units of 0.01 mm 





Deviations of tree no.— 





| Sum | Mean 
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Ficure 2.—Annual ring widths and tree trends of 12 dominant fir trees, 1864-1926. 
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TABLE 1.—Deviations of annual ring width from trends of individual trees, in 
units of 0.01 mm—Continued 


Deviations of tree no. 


Year rs eee Wace ie iene eet maemo | Sum | Mean 
1 2 3 4 ) 6 7 8 9 10 il 12 

1901 —21 12; -1 -13 40 9 | —13 36 | —11 —3 | -—49 | —27 | —41 -—3 
1902 0 31 30 8 42 2 58 38 30 | —2/ —10 11 | 288 24 
1903 10 41 71 68 65 24 68 70 61 19 10 | —10| 497 41 
1904 41 30 62 sy 57 37 89 53 | 72 40 20 18 | 608 51 
1905 —8 19 33 30 29 20 30 55 14 | —30 | —50 46 | 188 16 
1906 12 28 4 60 41 32 60 27 | 135 il 20 | 104| 584 49 
1907 43 |-112| —5 41 14 15 41 | —21 16 >| —21| -18| —5 0 
1908 23 37 25 31 26 38 31 22 87 23 19 40 | 402 33 
1909 14) —4| —24 |) —18 | —12 | —39 | —18 | —36 | —11 | -—16| -—1| —2 |—167 —14 
1910 24; 105 16 32 60 3 32 | —24 70 14 39 | 136 | 507 42 
1911 —6 4) -13 3 43 | —24 3 | —22| -19| —5| —41 | —66 |—143 —12 
1912 —26 14; —23 -17 15 | —51 | —17 | —40 12 | —24 | —32 22 |—167 —14 
1913 23 133 68 34 77 42 34 43| 123 37 68 | 130, 812 68 
1914 43 92 58 54 59 14 54 35 65 38 | 158 77 | 747 62 
1915 13 51 2 | -—6 12 27| -6| -3 16; 8] —42 45 | 143 | 12 
1916 42) 110 78 45 44 60! 45 29 87 19; 108} 122 | 789 | 66 
1917 22 10 | —5 6 133/ —5| —8/|-92| 10 17 10| —4 0 
1918 11 | —71 | —22 —55 | —11 | —25 | —55 | —26 0| —9| —23 | —23 |—309 —26 
1919 0; —12 8 -—44 -—29| -12| -44/ -4)} 1 12 7 | 4 |—113 -9 
1920 19 —43 | -33 -—4| —57 0| —4| —22| —58 | —18 57 | —18 |—181 —15 
1921 —22| -—54 -23 -—54 | -75| -—7| —54 | —20 | -—77 | —27 17 | —61 |—457 —38 
1922____ —33 -104 -—53 —54 | —82| —54 | —54  —57 | —95 | —36 | —64 | —94 |—780 —65 
1923 —44 | -75 -—34 -—24 | -—70 | —21 | —24 | —35 | —34| -—5| -—4 | —37 |-—407 —34 
1924. —36 | —56 | —14 6| —58| —8 6 | —13 | —43 26 | —54 | —30 |—274 —23 
1925... . —17 | -37| -15 -14 | —56 4|-14 29 18} 16] —54 | —53 |—193 —16 
1926... 21 | —27 | -—16 2 | —43 17 26 32 50 47 | —24 | —16 93 s 


MAGNITUDE OF ANNUAL GROWTH FLUCTUATIONS 
ANALYSES OF CURRENT ANNUAL GROWTH 


Dominant trees in an even-aged stand are free—and have prob- 
ably been free throughout most of their earlier life—from epochs of 
severe suppression by neighboring trees. As a consequence, progres- 
sive changes in their growth rates with respect to age (as distinct from 
changes with respect to time) should be subject to elimination by fit- 
ting trend curves, of first or second degree only, to the data of each 
tree. These should be sufficient to eliminate the effects of age on 
growth without danger of eliminating important changes in incre- 
ment due to climatic influence.* A straight line was therefore fitted 
to the values of trees nos. 2, 5, 6, 8, 9, 10, and 11, and a parabola to 
trees nos. 1, 3, 4, 7, and 12 (fig. 2). 

The residuals from these trends (table 1) were made the basis for a 
number of preliminary comparisons between the trees. Correlation 
coefficients between the annual growths of paired trees are as follows: 


Correlation 


coefficient 
Trees 1 and 2 0.422 
Trees 1 and 3 .633 
Trees 3 and 4 .510 
Trees 4 and 7 .310 
Trees 5 and 6 o .512 
Trees 7 and 8 34 .586 
Trees 9 and 10__-- 4 .594 
Trees 11 and 12 .608 


On the hypothesis that they express the relationship found between 
samples of independent variates, these correlation coefficients would 


‘ This statement is based on the assumption that if the 63-year climatic record can be expressed as a func- 
tion of time, the function is more complex than would be expressed by a curve of first or second degree only. 
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be distributed around zero with standard deviation of about 0.077. 
In the face of an average correlation coefficient above 0.5, however, 
such a hypothesis is untenable. On the other hand, with these mod- 
erate coefficients of correlation, there is a great deal of variation in 
annual growth between single trees unexplained by a system of causes 
outside the stand. If this unexplained variation is peculiar to the 
individual trees or to random error in measurement of the annual 
rings, it is subject to considerable reduction if the residuals from 
trend in each year for the 12 trees are averaged. 

As a step in this direction, correlation coefficients were calculated 
between average residuals of two sets of six trees each: 


Correlation 


coefficient 
Trees 1, 2, 3, 4, 5, 6, and 7, 8, 9, 10, 11, 12_- Seduce wa ee 
Trees 1, 3, 5, 7, 9, 11, and 2, 4, 6, 8, 10, 12__- - -. 945 


These high values demonstrate clearly that factors external to the 
stand have accounted for the bulk of the variation between the aver- 
ages of annual growths, based on as few as six trees each. The prob- 
ability that these factors are anything other than weather or climate 
seems extremely remote. 

The best estimate of the magnitude of variation in ring width 
that is due to fluctuations of climate will be arrived at from an analy- 
sis of variance (15) of the raw data of figure 2. From the classifica- 
tion in table 2, it is evident that a large proportion of the sum of 
squares of the deviations of each tree-year observation from the gen- 
eral mean of the annual ring widths of the 756 (12 times 63) observa- 
tions is due to differences between trees. Subtracting this quantity 
from the total, the remaining sum of squares is within the trees, and 
this in turn is divided into sums of squares (table 3) that are elimi- 
nated by the trends (age progression), by differences in years (effect 
of weather and climate), and by the residual, this last to be taken as 
the basis for determining random variation. Mean squares for these 
quantities are given in the last column of table 3. 

If the variation in the annual growth of these trees from year to 
year, after eliminating the effects peculiar to the individualities of the 
trees themselves and to changes in their ages, were entirely fortuitous, 
the mean square due to years would be expected to be practically 
identical with the mean square in residual variation; in fact, from the 
distribution of F as given by Snedecor (15, pp. 88-91) a ratio of the 
former to the latter higher than 1.25:1, when based upon the same 
number of degrees of freedom as in this case, is to be expected only 
about once in every 100 random samples. The observed ratio of 
about 15:1 is therefore taken as proof that the variation in growth 
from year to year is not the outcome of a constant system of chance 
causes. Again, it is seen to be due to factors external to the stand, 
= these can hardly be anything other than changes in weather and 
climate. 


TABLE 2.—-Division of the sum of squares into portions between and within trees 





a , ae Sum of Degrees of 
Source of variation squares freedom 

Between trees. ........ - 1, 328, 800 11 

| 2 : Ss atialciioie nee 1, 764, 480 744 


- a aesed kaeenae 3, 093, 280 755 
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TaBLE 3.—Division of the sum of squares within trees into portions due to trends, 
years, and remaining variation; and pertinent mean squares 

















Sum of Degrees of Mean 

Source of variation squares freedom square 
icons lsaeles tlie es dots iadbaasisiigi actu enisinds inrhasteioabeabalideamiimail 692, 336 117 40, 726 
NE ES LOE EL TE 620, 694 62 10, 011 
REREAD 2 451, 450 } 665 679 
WU oc cictsintctccutiinnnitinnmimttig ee nn een es 





1, 764, 480 | 744 








1 The 7 trends represented by straight lines account for 7 degrees of freedom; the 5 trends represented by 
parabolas account for 5X2=10 degrees of freedom; the sum of these is 17. 


The mean square between years is 12 times the variance of the 63 
annual ring widths (as measured from the tree trends) of the average 
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It contains one portion of variance due to the residual varia- 
The variance of annual ring width due to 


and its standard deviation is the square root of this quantity, or 27.8, 


in hund 


redths of millimeters. 


This is the magnitude sought. 


It may 


be taken as characteristic of the forest represented by these trees, 
because the growth of the entire forest is measured primarily in the 
growth of its dominant trees. 

Had the residual mean square of table 3 been zero, the variance 
between years would have been 


10,011 _ 


19 = 834, 


and it would have been concluded that, after eliminating the dif- 
ferences between trees and their age changes, fluctuations of climate 


account for the whole of the remaining variance. 
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The general mean of the annual ring widths during the 63-year 
period is 177.9 in hundredths of millimeters. 
due to annual fluctuations in climate is therefore 15.4 percent of this 
In 1 out of 20 years, accordingly, fluctuations in climate may 
be expected to be great enough to result in departures in diameter 


mean. 


growth of more than 30 percent of the longer time average. 


The standard deviation 
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variation in climate from year to year is itself the result of a mul- 
tiplicity of causes no one of which predominated, variation in growth 
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of periodic growth fluctuations of the stand represented by these 
trees which must be considered to be due to changes in climate is 


shown in table 4. 
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TABLE 4.—Coefficient of variation of periodic diameter growth 





(eeeey 
Deviation from 63-year 











Deviation from 63-year 
mean 
. . r | ‘ 7 3) | 
Length of period (years) Three Length of period (years) | | Three 
Standard times Standard times 
deviation standard || | deviation standard 
| deviation || deviation 
} | | 
Percent Percent } | Percent | Percent 
Be 15.4 46 || 15.- | 4.1 12 
5 6.9 | 21 || 20... | 3.4 | 10 
10 4.9 15 | } 


FLUCTUATIONS IN CURRENT ANNUAL GROWTH IN BASAL AREA AND VOLU ME 


The effect of fluctuations of climate on basal area or volume incre- 
ment may be deduced from two expressions: (1) The expression of 
basal area or volume in terms of diameter, and (2) the effect of fluctua- 
tions of climate on diameter increment. 

For the expression of both basal area and volume in terms of 
diameter we may take 

v=cd* 


in which d is tree diameter, c is a constant, and v is basal area if a=2, 
or volume if ais between 2 and 3. The value of a for volume deter- 
mination depends upon definition of volume, species, etc. 

If the annual diameter increment of a tree be denoted by i, and the 
standard deviation in 7 that is due to fluctuations in climate, ex- 
pressed as a ratio of i, be denoted by p (the value of 100p for these 
data being 15.4 as determined above), then 7 is subject to an error pi 
on account of fluctuations in climate. These definitions lead to the 
following determinations: 

(1) Without considering the fluctuations of climate: 


Ze 5 Diameter Basal area or volume 
Beginning of period TG a SR cd* 


End of period TRI MOC REE Gt: owcacstd 640 
Increment during period - . See, ee nae = 1 c(d+i)*—cd* 


(2) Considering the fluctuations of climate: 
J Diameter Basal area or volume 
enine GF MONGOE. coc casunaa d cd 


End of period , .------d+i+t pit c(d+i+ pi) 
Increment during period _- accents it pi c(d+i+ pi)*—cd* 
Increment fluctuation .____________- + pi ce(d+1t+ pi)*—c(d+1)¢ 
Increment fluctuation in percent of in- (d+it a. e 
d+t+pi)*—(d+i)° 
crement- ale ee eer +100p 100 (d+1)*—de 
Expanding the expression of increment fluctuation in percent of 
increment for basal area and volume by the binomial theorem, and 
retaining the first two terms only, in both numerator and denomina- 
tor, we obtain 





109) (@EO* tad +e (+ pi)+ -------------------- — (d+1)*]_ 
PUTS Gee ei —d* a 
a(d+7)*"!. (+ pt) 
of 2 apy 
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and this reduces to 
, “el 
100) 1 +5] - (+p) approximately. (1) 


If, for example, the average tree is 30 cm in diameter, and its rela- 
tive increment J is 0.10 in 10 years, the fluctuation to be expected in 


basal area increment in the next 10 years as a result of the caprice of 
climate is 
+ 100p(1+0.10)?-!= + 100p(1.10), 


that is, the expected relative fluctuation in basal area increment is 10 
percent greater than that of diameter increment. The expected 
relative fluctuation in volume increment, giving a the value 2.5 (2, p. 
279) would be + 100p(1+0.10)?*-'= + 100p (1.15), or about 15 percent 
greater than that of diameter increment. 

It is at once apparent that the relative fluctuation in both basal 
area and volume with respect to that of diameter depends upon the 


. U . . . . . . . 
ratio 7; that is, it varies directly with diameter increment of each 


tree—which, in turn, depends upon the annual growth of the tree and 
the length of period chosen—and inversely with the diameter of the 
tree itself. Hence, while real and positive, it is hardly sufficiently 
great to warrant its calculation for every tree, unless small timber of 
fast growth is being dealt with and longer periods than customary 
are used. The application of the relative fluctuation of the tree of 
average diameter (by basal area or by volume as the case may be) to 
the periodic growth of all the trees is sufficient approximation. 


ELIMINATION OF CLIMATIC FLUCTUATIONS IN MAKING COMPARISONS OF PERIODIC 
INCREMENT 


When periodic increments of a forest or sample plot have been 
determined from periodic inventories, comparisons between them are 
usually made for the purpose of estimating the effect on growth of a 
treatment, epidemic, or catastrophe, such as thinning, insects, or 
fire. The foregoing analyses, however, demonstrate that periodic 
changes in climate affect increment to a marked extent, independent 
of treatment or hazard. It may be necessary, therefore, to eliminate 
so far as possible the effect of climatic change in order to arrive at a 
valid estimate of the effect of treatment or other action. 

In the simple case of three inventories, or two periodic increments, 
it is not necessary to specify any particular method of increment 
calculation. Obviously the more accurate the measurements and the 
method, the more reliable the results. Only the adjustment of the 
periodic increments to the values to be expected if the effects of 
climate were the same during the two periods need be considered. 

The data upon which adjustment of this nature is to be based should 
be measurements of annual ring widths of 10 to 20 dominant trees of 
average size during a number of years equivalent to several periods 
during which the trees were relatively free from epochs of suppression 
by neighboring trees and from the influence of the treatment, epidemic, 
or catastrophe which is under investigation. If such trees cannot be 
spared from the sample plots or forest where periodic increments are 
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to be compared, they should, at least, be taken on the same species, 
same site quality, and in the immediate vicinity. Crosscuts along 
which measurements are made should, obviously, be above the in- 
fluence of butt swell. 

The data may then be subjected to analyses as discussed above. 
The effect of age may be eliminated by plotting the annual ring 
widths of each tree on time, and fitting a straight line or simple 
curve by freehand methods. It is not necessary to use the method of 
least squares, although with straight-line trends the work involved 
would not be great. 

Deviations from the trends may be classified as in table 1, according 
to tree and year, and the average (algebraic) deviation of the trees for 
each year computed. These may then be plotted on a graph as devia- 
tions from a base horizontal line representing the general mean annual 
ring width of all measurements. If the graduations of the vertical 
scale are doubled, the results, as shown in figure 3, express the average 





1880 1882 1884 1886 1888 1890 1892 1894 1896 1898 


TIME 
FIGURE 3.—Current annual diameter growth, corrected for age; periodic annual average for each of two 


periods, and general average diameter growth. 


of the annual diameter increments, independent of the individual 
trends. The graphs of figure 3 show: 


Mm 
Average diameter growth, first period. -...........------- 3. 85 
Average diameter growth, second period._.___....-_------ 3. 25 
Co ER Ed SE SR tae 3. 55 


Considering the diameter-growth fluctuation (due to climate) of 
the mean tree of the forest or sample plot as running parallel to that 
of these sample trees, it may be adjusted by the correction factors: 


oN) 
or 
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In order to make adjustment in basal area or volume increment, the 
diameter of the tree of average basal area or volume is first calculated. 
If the diameter of the tree of average basal area is 30.0 cm, or 300 mm, 


° . . 7 ° 
to a satisfactory degree of approximation, qd of equation (1) for a 10- 


vear period may be taken as 


Now we have first to calculate the factor p, i. e., the deviation of 7 
expressed as a ratio of mean 7; that is, 


p= +9. x a = +0.085, 


oo| HE 


The expected growth fluctuation then in basal area is, according to 
equation (1) 


100 (+0.085)(1+0.118)?-'= +9.5 percent. 


The basal area increment as obtained for the first period should 
therefore be reduced by 9.5 percent and the basal area as obtained for 
the second period should be increased by 9.5 percent in order to 
eliminate the influence of differences in climate in the two periods. 


INFLUENCE OF PRECIPITATION 
METHOD OF STUDY 


The methods used by Fisher (6) in studying the influence of rainfall 
on the yield of wheat in England formed the basis for the methods 
used in the present investigation. 

The problem is the calculation of the average influence of the rain- 
fall distribution during a period of several seasons, upon the width of 
the current annual ring of dominant trees. The length of the period 
should obviously be at least that of the growing season. It seems not 
unlikely, however, that the rainfall distribution during a number of 
months preceding the beginning of a season’s growth may have a 
measurable influence upon the width of the annual ring. There is a 
practical initial limit to the period, however, beyond which the in- 
fluence of rainfall upon the current season’s growth must become 
negligible. In the first investigations the initial limit was set at 
August 1 of the year preceding growth of the annual ring. The end- 
ing limit, which is taken as of August 31 of the year of growth, is 

rather definitely fixed by the time of year that the tree has completely 
formed the annual ring. The period taken is therefore 13 months. 

For purposes of the following discussion the term “regression func- 
tion” is defined as the function which expresses the influence of 1 ad- 
ditional centimeter of rainfall during a time interval upon that part 
of the width of the annual ring that is independent of the differences 
between trees and changes in age. It is denoted by y(t) wheret isa 

5 If to the data of the last column of table 1 were added the mean width of the annual rings of the 12 trees 
during the 63 years (177.9 in hundredths of millimeters), the results would be widths of annual ring adjusted 


for differences between trees and changes in age. These are the values referred to by the term ‘“‘width of 
annual ring’”’ or “‘growth.” 
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date within the 13-month period. Now there may be determined 
from the rainfall record (daily, weekly, monthly) during the 13-month 
period, another function, ¢(t), to express the distribution of rainfall 
during the period, from which the rainfall in any time interval t,—t, 
within the period would be given by 


te 
[ o(t)-dt. 


Because the influence of 1 cm of rainfall on the width of annual ring 
may be expressed by the regression function y(t) during the same 
time interval, the influence of the rainfall distribution on growth during 
this special time interval would simply be 


"te 


v(t) -p(t)-dt. 
If this integral is calculated over the entire 13-month period (the 


initial and final limit of which may be designated as — > and s) the 


whole effect of rainfall during the particular period upon radial growth 
may be written 
+2 


Growth= V=M- i y(t)-(t)-dt, 


-= 


M being a constant. This equation may be written for as many years 
as there are annual ring widths of the trees and 13-month rainfall «lis- 
tributions with which ring width is paired, the function y(t) repre- 
senting the average effect of 1 cm of rainfall to a unit of time (a day, 
a week, or a month) within the period, the unit of time being the same 
in all years. 
For the years 1, 2, 3, __n, we then have 
T 


V,;= -M+{ *¢ : t)-y(t)-dt 


ks 
“ 


a+; 
Vou + | », bo(t)-¥(t)-dt (2) 


2 


+3 
V,=M+ i 7 On(t)-w(t)-dt 
“3 


a a 


rene for ./ the arithmetic mean of ~ ring widths V;, V2, 


oe --- Vy, and letting 1. =(V,— M), m= (V2—M), ---- ---- 
v,= (Vv, _ ‘M), equations (2) become 
7 
Y= oi(t)-W(t)-dt 
’ ~<a 
we 
n=) r $2(t)-y(t)-dt (3) 
+f 


Uv, = p >, (t)-p(t)-dt 
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It is necessary, therefore, to’determine first the coefficients of the 
rainfall distribution, ¢(¢), for each of the 63 periods for which the pre- 
cipitation record is at hand, and then the coefficients of the regression 
function, y(t), such that the standard error of estimate of the n 
equations (3) is a minimum. 

The first step toward the solution of this problem is the choice of 
mathematical form of the functions ¢(f) and y(t). As to (ft) it is 
known that the distribution of rainfall during each of the periods of 
13 months is a waved curve. Knowing nothing else regarding it, the 
simplest function expressing such curves is the polynomial 


r1) (t) =) t+e,t+cet?+¢,f+ dl daca _-+e,t*. 


Concerning the regression function, y(t), it is known only that it must 
be a slowly changing function; a dry period of several days, for ex- 
ample, should influence diameter growth but little since the stems 
and the soil would contain enough water to balance the deficit for 
such a short period, It may, therefore, be expressed in the same 
form as ¢(t), that is, 


¥(t)=vot nitty? + --------- +t. 


Following the inquiries of Fisher (6), these rational polynomials 
may be written in the form of orthogonal polynomials. The develop- 
ment below will show not only that such transformations enable us 
to determine the necessary coefficients, but that other important 
advantages follow. Computations are very much reduced since 
orthogonal functions have been tabulated by different authors. 
Perhaps the most convenient tables are those of Lorenz (11), who 
presents an excellent discussion of the properties of orthogonal poly- 
nomials and the calculation of their coefficients. Only their appli- 
cation to the present problem will be explained here. 

A polynomial of rational function 


Y=Co ter! + Con? + - __--+¢,2 
may be written in the form of an orthogonal polynomial 


y = Ay 2 rD, ¢ — (yXo+ e a diene sod +a,X,, 


in which dp, a, @2, .--- a, are constants; and X,, which is a function 
of x that may be written X;,(r), is defined by the following determi- 
nant: 


rf 2z # Sa 

Tf ta! Sr Le 

r 7 ie = 9g ~ 4 | 
X,=b4Ze Ze Te... ......2e" | 
Ye? De drt. Dat? | 


¥,i-1ly ,t ¥ ptt! 
at ae” seseen 
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Properties which make these functions orthogonal are the following, 
in which i¥k: 


- 
>X,=0 or | *x-ar=0 
/T 
wr 
SX,-X,=0 or | °X,Xed2=0 4) 
T 
l SX? =] or 4 XY? -dr= 
X= X?,.dz= 
N or 


The constant 5, in the determinant is so chosen that >X?,;=N, 
where N is the number of equidistantly spaced observations of z. 


From the coefficients a, @;, d2, ---- a, of an orthogonal polyno- 
mial, the corresponding coefficients ¢», ¢, C2, ---- ¢,, of the poly- 


nomial of rational functions may be calculated, if needed. The latter, 
however, are not necessary to the present analyses. 

The application of orthogonal polynomials in the present investiga- 
tion consists of the following steps: 

(1) The expression of the rainfall distribution within each time 
interval of each of the 13-month periods by the orthogonal polynomial 





o(t)- = Ag+ a,X,+a,X.+ oaocesecee + Ope ¢ (5) 


for which the determination of do, a, a2, -.-. a, will be illustrated 
in the next section. 

(2) The calculation of the regression function, which is to be 
expressed by another orthogonal polynomial 


v(t) =agt+aX1+aX2+--------- + aX; (6) 
for which the method of determining ap, a, a2, ---- a,, Will next be 
discussed. 

According to equations (5) and (6) we have 
vy +? 
Vi | r v(t) ° o(t) -dt= f r (dg +a,X,+ ee a,Xx) (a+ aX, 


2 


whence, on account of properties expressed by equations (4), in which, 
when i+k, 


7 
f 
sat 


as . 
2 , r al r q 
X,:X,-dt=0 and X?, -dt=T, 
rT Bs 
2 


it follows that 


0,= T (argo + 040, + Op0g+ --------- + aay) +C (7) 


where C is the constant of integration.° 


6 Footnote on opposite page. 
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After we have obtained the rainfall coefficients dp, a,,-..____.- x, 
for each of the sixty-three 13-month periods, we may set up 63 observa- 
tion equations of (7), the unknowns, ap, a,--------- a,, of which are 
the coefficients of the regression function y(t). And these may be 
calculated by the method of least squares. The multiple correlation 
coefficient obtained at the same time would then afford a measure of 
relationship between growth and rainfall distribution. Thereafter, 
upon multiplying the a, coefficients by the corresponding functions 
X, of equation (6) we obtain the regression function y(t) which gives 
the effect of 1 additional centimeter of rainfall per month upon the 
width of the annual ring. 


RAINFALL COEFFICIENTS AND THEIR PROGRESSIVE CHANGES 


The rainfall distribution coefficients for each of the periods might 
be calculated from the rainfall record according to whatever time 
interval—daily, weekly, monthly—the record is preserved. Little, 
though, if anything, would be gained by using an interval which from 
phenological observations is less than that between calendar dates 
which mark the earliest activity of the buds in spring. Roughly, this 
interval is 2 weeks. The record of rainfall at hand, however, is accord- 
ing to successive months only. The span of a month seems long, and 
the rainfall distribution coefficients may, in consequence, not be as 
accurate as might be wished. 

The first polynomials were fitted to the monthly precipitation 
values for a 13-month period beginning with August 1 of the year 
preceding growth and ending with August 31 of the year of growth in 
each of the 63 years. These were of the sixth degree. The calcula- 
tion of do, a;,-...-----ds, for the annual-ring year 1924, is illustrated 
in table 5. The symmetrical arrangement in equal time intervals, t, 
measured from the median interval, is essential if the tables of orthog- 
onal functions prepared by Lorenz (1/1) are used. The form of table 
5 (known as the scheme of Chetverikoff) is that suggested by Lorenz 
when N is odd. The coefficients do, a, d2,---.----- d, are computed 
through the use of cumulative sums, combined with certain parame- 
ters which have been computed by Lorenz (11), by whom the tech- 
nique of these computations is fully described and the necessary tables 
given. 

The coefficients dp, a, d2,--------- a, are those of the orthogonal 
polynomial equation (5), according to the method of least squares. 
An immediate advantage of the orthogonal polynomial form is that 
coefficients up to a given degree are not affected if coefficients of higher 
degree are added. By simply adding the term a;X; to the polyno- 
mial of second degree we obtain the polynomial of third degree, or by 
dropping the a,X, term, the remainder is the polynomial of first degree. 

6 Certain remarks concerning C and T may prove helpful: 

The constant of integration must be C=—T' (aofo+-and;+--------- +ardx) where do, 4;,..........-de are 
the means of the corresponding rainfall coefficients of the n years considered. These mean coefficients 
characterize the average rainfall distribution of the n years, and consequently the deviation of growth 
6; for average rainfall conditions should be zero, because the sum of the n equations (3)=0. It follows 


necessarily that 
5;= T (cxodio-andy+ --...-- --tariz)+C=0, 


if Cis the quantity given above. 

Since 7’ is the entire time period under consideration (in this case, 13 months), and since the rainfall 
coefficients, ao, ai, @2, -.-....-.- a», are calculated on the basis of rainfall data obtained for a certain num- 
ber, N, of time intervals AT (using the functions X; corresponding to this number of observations N), 7’ 
must be N-AT’.. In this discussion AT=1 month, hence T= N-AT=13. 
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TaBLe 5.—Scheme of Chetverikoff for calculation of rainfall coefficients, annual 
ring nd bind . 











Month t w 
ee we = Sa _ — 

1923 | 
August _- —6 | 51 51 51 | 51 51 | 51 51 51 
September..|—5 | 80 131 182 233 284 335 386 437 
October —4 |227 358 | 540 773 1, 057 | 1, 392 | 1,778 | 2,215=s8-5 
November__|—3 /141 499 | 1, 039 | 1,812 2,869 | 4, 261=s-; | 6,039=s-; 
December. _|—2 |172 671 | 1,710 | 3, 522=s8-2 | 6,391 = 8-5 ‘ 

1924 
January —1 | 38 709 = 8. o | 2,419=s8-; 
February 0 | 26 26= uy 
March 1 | 55 717 =a 2,596=3a, | 6, 797=82 
April 2 |138 662 1, 879 | 4, 201 8, 130=8, | 14, 258=8, 
May 3 |153 524 | 1, 217 | 2, 322 3, 929 6, 128 | 9, 009=s5 | 12, 662=s¢ 
June 4 |139 371 693 1, 105 1, 607 | 2, 199 | 2, 881 3, 653 
July 5 |142 232 322 412 502 | 92 | 682 772 
August 6 | 90 90 90 90 90 90 

8-0-+-U'9+-80= So |81—8-1 = S; |82+8-2= Se |83—8-3= Sy |84-+-8-1= Sy \85—8-5= Ss |85-+8-0= Se 

Ss 1, 452 177 10, 319 1, 739 | 18, 519 2,970 14, 877 


a ; ao= 111,692 | a;=3. 639 | a2=0.825 | a:=11.956 |a,= —42.691 as= —0.058 |ag= 24. 155 


1 w=rainfall in millimeters. N=13 


Upon combining the coefficients of table 5 with the corresponding 
saan functions, as tabled in Lorenz’s paper (//), the values of 
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FiauRE 4.—Monthly precipitation from August 1923 to August 1924 compared with the polynomial curves 
of fifth and sixth degree. 





the precipitation curve are at the points corresponding to the months 
for the annual-ring year 1924. In figure 4 the curves of both the 
fifth and sixth degrees are presented. As may be checked in table 5, 
the first coefficient, a, is simply the arithmetic mean of the monthly 
rainfall records. The second coefficient, a,, leads to the average 
change of rainfall per month during the period, the third coefficient, az, 
measures the parabolic term, and so on as more complex features 
of the precipitation distribution are taken into account. 

The coefficients for the 13-month periods are listed in table 6, and 
they are presented graphically, by 7-year moving averages, in figure 5. 
Aside from whatever bearing these coefficients may have on tree 
growth, they are of interest in themselves in that the moving averages 
of figure 5 suggest systematic changes in their values during the 63 
years. In order to test whether these are significant, a polynomial 
of sixth degree, of the same form as equation (5), was fitted to the 
63 observed values of each coefficient of table 6. These, in turn, 
are listed in table 7, the first line of which—the a’, coefficients—are 
the averages of a, from the corresponding column of table 6. 
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FIGURE 5.—Polynomial coefficients according to 7-year moving averages. 
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TABLE 6.—Coeffcients of orthogonal polynomials which express teen trends 
(in millimeters) by monthly intervals from August of a preceding year to August 
inclusive of a current year, for the current years 1864-1926 





Coefficients for— 





Current year — — ee 








ao a a2 a3 a, as ae 

Ist 80 -3 10 —16 —17 -1 14 
1865 73 22 | 14 12 14 0 -4 
1866 100 24 27 -1 39 12 21 
1867 118 —13 9 —28 24 —28 16 
1868 71 -3 16 4) —5 4 —13 
1869 &S 3 13 12 | —N -1 19 
1870 74 | —2 38 12 | 3 —5 s 
1871 94 | —ll 22 0 | —19 —5 12 
1872 M4 28 17 —18 | —12 12 21 
1873_. x 100 ll 12 —17 | —-14 —15 20 
1874_.... 70 9 22 —11 —-2 2 6 
1875 . : 81 Ss 5 8 —12 —30 14 
eae . 111 6 —15 —15 | —12 10 —2 
1877 89 14 1 —8 2 3 —15 
1878 110 27 14 —2 | —ll1 —8 15 
1879 101 7 21 —7 | 2 —22 —5 
1880 80 2 27 0 8 7 36) 
1881 a 94 —5 31 12 4 20 12 
1882 s 91 —2 61 —10 —2 7 -4 
1883___. a 114 —30 37 4 —34 —20 —6 
1884 73 6 24 14 —12 —2 3 
1885 65 —1 12 —22 0 -7 ll 
1886 " 95 20 31 7 —30 1 —9 
1887 on 81 —1 24 1 -3 —14 18 
1888 80 18 19 —8 5 -9 -Y 
1889 sO 1 31 -—7 —20 1 —18 
1890 81 20 41 19 2 12 20 
1891 _ 81 —2 46 —24 -3 —23 12 
1892 ‘ s4 4 13 2 —10 —15 —2 
1893 69 —13 14 —Ss —20 2 —19 
1894 74 20 i) —11 —25 10 —3 
1895 73 0 7 0 -4 1] 0 
1896 90 25 16 7 0 -—3 17 
1897 102 —2 44 2 13 29 —10 
1898 ot | —5 26 —5l —4 -4 —5§ 
1899 69 6 6 1 —13 2 16 
1900 80 6 12 9 1 —2 —22 
1901 92 19 14 —10 1 —*s 2 
1902 SS 3 15 —12 3 5 —10 
1903 74 11 28 5 -—2 -3 —3 
1904 = 90 —5 7 —11 6 4 10 
1905 92 31 31 3 14 13 4 
1906 at 4 —18 22 -—31 3 —16 ; 
1907 ne 75 1] y 3 —2 —14 —6 
1908 0: SS 13 4 -—6 —4 -9 —13 
1909 M 20 24 —5 —17 —* —18 
1910 130 26 10 4 —10 —39 —11 
1911 M4 —24 10 —6 -l -9 22 
1912 108 36 19 12 —5 | —2 4 
1913 a 97 —22 25 —s 3 -9 13 
1914__ nee 122 27 14 2 0 1 —8 
1915 103 2 17 —18 13 —14 —7 
1916 114 13 11 3 —8 —22 —8 
1917 113 16 30 ot) -4 —4 8 
1918 93 -—15 30 —18 —17 Rk 10 
1919 97 —14 —12 0 —7 -3 —15 
1920 M —2 —15 33 —ll —32 19 
1921 60 7 20 —6 cat | Q 2 
1922 108 14 -1 —12 14 —2 —s 
1923 93 —23 —6 ~a4 -1 3 ll 
1924... 112 4 1 12 —43 0 24 
1925_. 79 17 26 0 —6 9 —4 
1926_. 100 2 0 —21 —10 —32 2 

ES Sea eer re 90. 1 5.5 16.7 —3.5 —4.9 —4.1 0 
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TABLE 7.—Coefficients of orthogonal polynomial equations which express trends of 
the a; congeners of table 6, aren the 63 years 1864-1926 








Coefficients of polynomial equation for the following coefficients from 
table 6 


Coefficient 


a’o.- . ane 90. 100 5.500} 16.700) —3.500 
3. 142 —.605 | —3. 586 . 322 

2. 518 —.743 | —3.983 . 743 

ae 7 —.742| —1. 250 1.101 | . 334 
A i's Gases —4. 503 —. 040 .572 | —1.455 
.... ; ‘ —. 043 1.626 | —2.070| —1. 936 | 
4 a's. : 2. 765 1. 507 -663 | —1.129| —.998 | —2. 
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The test is carried out through the analysis of variance of each a; 

coefficient of table 6. In general, 
=(y—g=2(y—Y)?+2(¥—g)? 

in which y is a dependent variate, 7 is the arithmetic average of these 
variates, and Y is the trend value of y as estimated by the method of 
least squares. In the analysis of variance of a as an illustration, the 
left-hand member of the above equation is the sum of the squares of 
the 63 deviations of a, from their mean, as computed from the values 
of the a) column of table 6. It may be written 








Y(ay>—a’y)*= Za 2,— 63a’, 
in which a’, is the arithmetic average of the 63 observations of dp. 
The term corresponding to =(Y—g)? may, on account of the orthog- 
onality of equation (5), be simplified to 
63(a’*,+a”",+a”",+a",+a";,+a"s), 

obviating the need of calculating the trend values of a) (11, pp. 44-44). 
The sum of squares independent of the trend [=(y—Y)’] 1s then the 
difference between the total sum of squares [=(y—¥)*] and the sum of 
squares due to the trend [[(Y—9)’]. Calculated thus, the sums of 
squares for the analyses of variance of do, a;, d2,-_--.---d,, and their 
+ degrees of freedom, are given in the second, third, and fourth columns 
of table 8; mean squares and the observed values of Snedecor’s (15) 
ratio F are in the last three columns. 


TABLE 8. nang of variance of rainfall eonpiarente 
Z 


Sum of squares and degrees of 
freedom in— 
Deviations 
from trend, 
56 degrees 

| 





Mean square in— 
Coefficient 


Deviations 
from trend | 


Total, 62 
degrees 


Trend, 6 


Tren¢ 
degrees end 


469. 35 204. 08 
77. 70 221. 06 
366. 45 208. 01 
. OE 188. 55 | 
270. 90 160. 09 
162. 7f 172. 35 
194. 25 165. 38 


ao 5 2, 816. 
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at 5-percent level of significance. 


27 
15 at 5-percent level of significance. 


1 F=2., 
a F=3. 





From the distribution of F in random samples when based on 
6 and 56 degrees of freedom as in this investigation, the value of F at 
the 5-percent level of significance is 2.27, provided the greater mean 
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square is based on 6 degrees of freedom, and 3.15 if based on 56 degrees 
of freedom (15, table 25). The observed F for ap is the only one that 
exceeds expectation at this level of significance. With this single 
exception, then, it appears that there is no real evidence of systematic 
changes of the tainfall coefficients during the 63 years with which this 
investigation is concerned. In successive years a given coefficient— 
other than ay—seemingly takes random values from a distribution 
whose mean and standard deviation are given in table 9. The mean 
and standard deviation of a, are also listed in the table. The mean of 
each coefficient is taken from the first line of table 7 and the standard 


deviation is the square root of \% times the sum of squares of the 
totals in table 8 


Taste 9.—Means and standard deviations of the rainfall coefficients, 13-month 


periods 
a a a) a a a as 
Mean 90. 100 5.500 | 16.700, —3.500| —4.900 —4.100 3. 000 
Standard deviation 15.157 | 14.304 | 14.945 | 13.354 13.069 13. 093 12. 968 


The mean values of these coefficients differ significantly from zero 
owing to the variations of monthly rainfall within the period. They 
represent the average sequence of rainfall, presented graphically in 
figure 6 in comparison with the actual monthly averages. 
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vy Y 
Preceding Year Current Year 
FiGune 6.--Average monthly precipitation from August to August, inclusive, and its representation by the 
average of the polynomial coefficients to degree 6. 


INVESTIGATIONS PRELIMINARY TO THE REGRESSION OF GROWTH ON RAINFALL 


It has often been observed that two variables, known to be intrin- 
sically independent of each other, are apparently correlated as a result 
of parallel progression of their observed values in time. Such rela- 
tionships are known as spurious, or nonsense, correlations (2, pp. 
287-289). An estimate of the true correlation between variables 
which are undergoing progressive changes in time may be calculated 
from the residuals from secular trends, provided the trend curves 
successfully eliminate secular changes. It must be emphasized that 
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the true form of secular trends is not usually known. While a trend 
curve of polynomial form may be fitted to time series with any 
desired degree of approximation by simply increasing the degree of 
polynomial, yet the decision as to the sufficiency of a given degree is 
not universally and entirely based on objective analyses, even though 
complete use be made of statistical checks. The decision is, rather, 
a common sense one, in that intimate acquaintance with the data 
under investigation is a most important consideration in the judgment 
as to whether a given trend is reasonable. 

Attention is called to these difficulties because the following dis- 
cussion goes beyond the results of mechanical application of certain 
well-known statistical formulae into a realm that may not be entirely 
free from the idiosyncrasies of individual judgments. 

It was pointed out in connection with table 8 that the trend of the 
rainfall coefficients ag—the average monthly rainfall from August of 
a preceding year to August inclusive of a current year— could hardly 
have been the outcome of chance fluctuations of annual rainfall dur- 
ing the 63 years. Consequently the question arises as to whether the 
best estimate of the effect of 1 additional centimeter of rainfall per 
month on tree growth is to be obtained from the observed rainfall 
coefficients a), or from the residuals of a) from its secular trend. The 
tree data also, while they are already averages of residuals from 
separate trends designed to eliminate progressive changes peculiar to 
ach of the 12 trees, may still be subject to collective progression, 
common to all the trees in the stand; such, for instance, as may be 
due to changes in soil conditions over the 63-year period. 

In order to eliminate progressive changes both in the tree-ring widths 
and in the rainfall coefficients a,, trend curves as expressed by orthog- 
onal polynomials of degree 1 to 8 were fitted to the annual observations 
of both variables, rendering unnecessary the tedious computations of 
the 88 sets of residuals of each variable upon which the direct calcu- 
lation of the correlation coefficients would be based. These correla- 
tion coefficients may be expressed 


. , 
LYK’ ks 


Por. : 72 
VU XY ev 


Vreks 


(8) 


in which y is a residual of an average annual ring width from trend of 





degree k [k=0, 1, 2, - - - - - 8] as indicated by its subscript, and y’ 
is a residual of average monthly rainfall during the 13-month period 
from trend of degree k’ [k’=0, 1, 2, - - - - - 8] also indicated by its 


subscript. But following Lorenz (11, pp. 49-50) the calculation of 
the correlation coefficients may be greatly shortened, for equation (8) 
may be written 


TT 
~Yoy’o— Nn (agg +a", + po ge est +a,0’,) (9 
\ [Syo’ —n (ay? +4," + ----- +-a,”)|[Zy’o?—n(a'9? +0’ 7+ -----+4’,/7)] ~ 
in which dp, a, - - -- - a, are the coefficients of the orthogonal 
polynomial of degree k, fitted to 63 average annual ring widths y, and 
lp’, Q;', ----- a’,, are the coefficients of the orthogonal polynomial 
of degree k’, fitted to the 63 periods of 13-monthly rainfall averages. 
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These coefficients have already been calculated as steps in their re- 
spective polynomials. The only new computation is the single sum 
of the products (y, y’,) common to all the correlation coefficients dis- 
cussed below. The numerator of equation (9) for a correlation coeffi- 
cient between residuals from trends of unequal degree is carried out to 
the higher of the two degrees, while the terms in the denominator are 
taken to the degree of the trend with which the residuals are identified. 
Writing 7 as an illustration, 


LYoYo’ —N (Ago’ + aA,’ + aod’) 


v [2 yo? —nag?][Z yo"? —m (ap"?+-a,"?+-a2”*) | 





ro2= 








Using equation (9) the correlation coefficients are as listed in table 
10. The first column and first line of correlation .coefficients are of 
particular interest, the first column containing those between mean 
rainfall (i. e., residuals from trend of degree 0) and the average residu- 
als of annual ring width of the 12 trees from trends of degree 0 to 8. 
The highest correlation coefficient in this column is 73, between ring 
width residuals from trend of degree 3 to rainfall, the value of which is 
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FicuRE 7.—Annual ring width data and their secular trend as represented by the polynomial of degree 3. 


0.367; between residuals of ring width from trends of higher degree 
than 3 and rainfall, the correlation coefficients fall off to 0.243 for rgp. 
Evidently trend curves of degree higher than 3 cannot be regarded as 
eliminating only the secular changes in tree growth; they also elimi- 
nate part of the annual growth fluctuation that is dependent on rain- 
fall, as indicated by the generally declining values of the coefficients 
from rg to 7. On the other hand, if one may judge by the increasing 
correlation coefficients between rainfall and annual ring width residu- 
als from trends up to degree 3, these lower trends eliminate systematic 
changes in ring width that have nothing to do with rainfall; hence the 
annual ring width residuals from trend of degree 3 are probably the 
best values for use in the regression function. The annual ring width 
data and the third-degree trend fitted thereto are presented in figure 7. 
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There are analogous conditions in the first line of table 10. The 
correlation coefficient ro.—between annual ring width and rainfall 
residuals from trend of degree 2—is the highest with a value of 0.374, 
dropping gradually to 0.218 in rg. Through the trends of degree 
higher than 2, rather short fluctuations in annual rainfall are elimi- 
nated, and these evidently affect annual growth. In the light of all 
the correlation coefficients in the table, it seems at first that the best 
estimate of the effect of an added centimeter of rainfall on tree growth 
should be obtained from the correlation between the residuals of an- 
nual ring width from trend of degree 3 and the residuals of the dp 
coefficients from trend of degree 2. However, since the increase in 
the correlation coefficient ro [=0.367] over ro [=0.348] is considerably 
greater than of 3; [=0.380] over ro. [=0.374], but little would be gained 
through the use of residuals of the rainfall values a) from trend, rather 
than the a values themselves. For it seems likely that the trend of 
degree 2 through the rainfall coefficients a, eliminates part of the rain- 
fall which has influenced tree growth. Furthermore, spurious corre- 
lation is avoided if the secular trend of but one of the two variables is 
eliminated; hence the original rainfall values are used for the inde- 
pendent variable dp. 


TaBLE 10.—Correlation coefficients between residuals of annual ring width and of 
annual rainfall, taken from trends of degree 0 to 8 





, : : Rainfall variable:residuals from trend of degree indicated 
Annual ring width varia- | 
ble:residuals from trend |- —— _ 





of degree indicated 














0 1 2 3 4 j 6 7 8 
| | | 

0 niles ---------| 0.348 | 0.355 | 0.374 | 0.357 0.299) 0.27! 0.234 | 0.219 0. 218 
a . , : . 347 . 355 . 374 . 357 . 299 . 275 . 234 . 219 218 
a . 362 . 370 . 375 . 358 . 300 .301 | .235 . 219 219 
eee a R 367 | 375 . 380 -381 | .319 . 320 + 250 . 233 
4 302 | . 314 314 | 330 . 331 . 259 . 241 241 

311 | 318 - 323 . 323 | 340 . 340 - 266 . 248 248 
| ee ae ee meh . 269 . 270 284; .284/ .200 . 270 . 270 
_ ae : 243 | .249 . 252 . 253 . 266 266) .271 274 | . 274 
a e ‘“ . 243 248 | .252 | 252 | 265 | =. 265 271 274 . 274 





REGRESSION OF GROWTH ON RAINFALL AT DIFFERENT SEASONS 


As an outcome of the analyses just discussed, two regression equa- 
tions were calculated, each expressing the average influence of the 
rainfall coefficients a», a,,-.-------- a; for the 13-month periods on 
annual ring width. The a, coefficient was not used as an independent 
variable because it seemed likely that the others should suffice to 
describe the influence of precipitation adequately. The observa- 
tions upon which these regression equations are based differ only in 
the values of annual ring width. In the first, the annual ring width 
residuals from trend of degree zero—the values of table 1, denoted by 
the symbol v—\were taken; in the second, the residuals were taken 
from trend of degree 3—that is, the deviation of the values of table 1, 
last column, from the trend curve of figure 7, denoted by vo. 

Observation equations of the form of equation (7) in which 


vp = T (ado ay0,+ --.------- + ass) 


were written down for each of the 63 years, the vo values for which 
were taken from table 1 and the values of the a; rainfall coefficients 
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from table 6 for corresponding years. These equations lead to the 
normal equations, which when solved simultaneously, render 13— 
since 7’=13—times the coefficients a, of the regression function. 
Then upon putting the a, coefficients in equation (6) in which 


V(t) =a + .X1+ 02X24 ---------- +asX5, 


where the functions X, are taken from Lorenz’s tables (11), we have, 
finally, the expression of the regression function y(t). This is pre- 
sented graphically in figure 8. 

In like manner, the residuals of annual ring width from trend of 
degree 3 lead to the regression function y(t) of ce; and this is also 
presented in figure 8. 
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FiGURE 8.—Average effect of 1 additional centimeter of precipitation per month on width of annual ring 
(13-month period). The curve marked e® is based on the annual ring width data of table 1 without 
further adjustment, while that marked v@ is based on deviations of the same data from the trend curve of 
figure 7. 


The two regression curves are quite similar. With the exception 
of March and April, an extra centimeter of precipitation per month is 
beneficial to the growing tree, although with varying efficiency in the 
different seasons of the period. The most striking feature, common 
to both curves of figure 8, is that the efficiency of an added centimeter 
of precipitation during the months of the year preceding growth is 
apparently as great as that during the season of growth itself, a higher- 
than-average precipitation during the late summer of the preceding 
year contributing most to the current season’s growth. During win- 
ter the efficiency of extra precipitation falls off, perhaps because it is 
largely in the form of snow, and scarcely available, if at all, before the 
spring thaw, for in early spring the curves reach a minimum. They 
then pick up steadily through May and June, reaching a maximum in 
July. Thereafter they decline with the completion of the season’s 
growth. 

The negligible influence of extra precipitation during August on 
the current season’s growth, as contrasted with its relatively great 
influence on next year’s growth, suggests that August precipitation, 
and perhaps that of the earlier summer as well, is especially stimulating 
to the production of food reserves needed for next year’s growth. 
But as the regression curves of figure 8 hang loose-ended relatively 
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high above zero in August of the preceding year, the 13-month period 
apparently does not describe the influence of precipitation on width 
of annual ring with sufficient completeness. 

It was for this reason that the rainfall coefficients and the regression 
function were finally calculated for the 15-month period, thus including 
the June and July precipitation of the preceding year. Furthermore 
the a, rainfall coefficients for the 15-month periods were used as a 
seventh independent variable. Two regression curves—one for v 
and one for v“’—were accordingly calculated, and these are presented 
in figure 9. 
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FIGURE 9.—Average effect of 1 additional centimeter of precipitation per month on width of annual ring 
(15-month period). The curve marked vo is based on the annual ring width data of table 1 without 
—eOO while that marked 0{) is based on deviations of the same data from the trend curve 

This time, as anticipated, the regression curves decline toward zero 
on the left-hand end. Otherwise the general trends are much like 
those of the 13-month period. The area under the curves—particu- 
larly under v®, which is based on annual-ring observations from which 
progressive changes not affected by precipitation have been elimi- 
nated—represents the relative efficiency of 1 additional centimeter of 
precipitation per month. The greater efficiency of additional precip- 
itation from June to December of the year preceding the formation of 
the annual ring over that of the 8-month period of the current year is 
again apparent. 

But perhaps the most interesting feature is brought out in figure 
10, in which the curve of v for the 15-month period is broken into 
two parts and set along the time scale of a year from approximately 
the initiation of one season’s ring to that of the next. Curve a repre- 
sents the influence of an additional centimeter of precipitation per 
month upon the growth of the current season’s ring, while curve 6 rep- 
resents its influence on next year’s ring. The new point of interest is 
that the influence of extra precipitation on next year’s growth becomes 
really pronounced at about the same time (early July) that its influ- 
ence on this year’s growth passes its maximum, and reaches its own 
maximum, in turn, at the end of this year’s growing season. During 
the month of July, primarily, additional precipitation contributes to 
the current season’s active growth, and at the same time apparently 
stimulates the building up of reserves in preparation for the spring 
get-away. From August on, the preparation for next season continues, 
although at a lessened rate. 
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No curve, however, may be regarded as representing exactly the 
true influence of the precipitation distribution on growth. As in all 
correlation analyses, the curves of figures 8 and 9 hold within certain 
limits of error only. Graphically the limits of acceptable error might 
be portrayed by a painted strip covering the two curves of each figure 
the trend of which should probably be guided mostly by the curve 
of »®. 

The sum of squares of the dependent variable is reduced in the 


proportion a Rey Where R is the multiple correlation coefficient. 


But observations during 63 years only give 62 degrees of freedom in 
e to start with; while in » there are but 59 because 3 degrees of 
freedom have been eliminated by the trend of third degree, the re- 
siduals from which are »™ by definition. At most, therefore, the 
samples are not large, particularly if the correlation coefficients are 
moderate. 
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FIGURE 10.—Average effect of 1 additional centimeter of precipitation per month during 1 year on (a) the 
current annual ring and (6) next year’s annual ring. 





Table 11 contains the analyses of variance of the data based on the 
regression equations of the 13-month periods. The sum of squares of 
both » and v® are divided into parts due to, and independent of, their 
respective regression equations, the former accounting for 6 degrees of 
freedom in each case. The correlation coefficient between actual and 
: et @) 10,168 
calculated values of »™ is fis 
correlation between growth and the rainfall distribution. The ratio 


ae 
oe =2.28. From Snedecor’s table (15) the value of F at 
‘ 








=0.443, indicating only moderate 


F of »™ is 2 


the 5-percent level in the case of independent variates based on the 
same number of degrees of freedom as are available here, is 2.27. 
The regression is therefore adjudged significant. 

The correlation coefficient between actual and calculated values of 


» ;, /10,272 ; ae Se . 
v® is *—- ==().476 and the observed ratio F is ——-=2.59. For 
? J 45,342 i 662, °° 

the same number of degrees of freedom the value of F at the 5-percent 
level in the case of independent variates is 2.28. This regression, 


consequently, is also judged significant. 
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TABLE 11.—Analyses of variance based upon the regressions according to 13-month 
and 15-month periods 


13-MONTH PERIODS 





} 
| Sum of |Degrees| 














Degrees 
S sass Mean _ Tr Sum of 1 Mean 
Source of variation : Source of variation | : 
SqUATES |freedom| squares || | squares |r .aqom| Sauares 
| 
aan een ielbeeeees ‘yed ‘mie: caer) \ ene 
Analysis of o® | Analysis of v0: 
Regression...........| 10, 168 | 6 1, 695 Regression...........| 10, 272 6; 1,712 
Deviation from re- | Deviation from re- 
gression..........| 41, 556 | 56 742 gression..........| 35, 070 53 662 
0 | 51, 724 62 ineitins ee er 
15-MONTH PERIODS 
Analysis of 0: | Analysis of 0(): 
Rogression...........| 8, 463 7 1, 209 Regression._..._- ----| 9,044 | 7 1, 292 
Deviation from re- Deviation from re- 
GUNTER cnccccss | 43, 261 | 55 787 REE 36, 298 | 52 698 
_——EE ——- | 





i ot | 51, 724 | 62 |... | ie 45, 342 | 59 





Table 11 contains also the analyses of variance based on the 
15-month periods. While this greater period length and the inclusion 
of the additional independent variable—the a, coefficient—are 
believed to supply a more complete description of the influence of 
rainfall on annual-ring width, a certain loss of significance results, if 
the latter is placed at the 5-percent level; for the value of F at this 
level, taken from Snedecor’s table, is 2.19 for the 7 and 55 degrees of 
1,209 
787 
sponding values of F for the 7 and 52 degrees of freedom upon which 
the analysis of v is based are 2.20 and 1.85. Several reasons, other 
than the fact that the size of sample is not large, may be suggested 
for the lack of positive proof of the calculated influence of the rainfall 
distribution during the 15-month periods on the timber represented 
by these data. Among them are the following: (1) The rainfall 
record according to intervals of 1 month may cover up part of the 
correlation that weekly or fortnightly records might bring out. (2) 
The trees from which the annual ring width measurements were 
taken grew at a distance of about 8 km from the weather station, so 
that the rainfall record used may not have been precisely that to 
which the trees were subjected. (3) In a region where the average 
monthly precipitation is 9 cm, rainfall distribution throughout the 
year is probably far from a limiting factor in timber growth. 





freedom for v, whereas the observed F is =1.54. The corre- 


SUMMARY 


The annual ring widths from 1864 to 1926 of 12 dominant Abies alba 
sample trees taken from a 120-year-old stand of fir-spruce-beech, and 
rainfall records by monthly intervals for this period, were available 
for a study of the effect of climate on the magnitude of timber-growth 
fluctuation and an investigation of the influence of rainfall distribu- 
tion on tree growth. 

After plotting the annual ring width on the year of growth for each 
tree, trend curves of growth, of first or second degree only, were fitted 
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to the data of each tree in order to eliminate the effect of age on width 
of annual ring. 

An average correlation coefficient of about 0.5 was found between 
residuals from trend of single trees, and of about 0.9 between average 
residuals of six trees from their trends, proving that an important 
source of variation of annual ring widths lies outside the forest. 
This source is taken to be the annual fluctuation of climate. The 
magnitude of variation of annual ring width due to fluctuation of 
climate as calculated through the analysis of variance is a standard 
deviation that is 15.4 percent of the average ring width of the 12 
trees during the 63 years. 

The effect of fluctuation in climate on basal area and volume 
growth are calculated, and a method of eliminating the effect of 
climate in comparing the effect of treatment, epidemic, or catastrophe 
on periodic growth is illustrated. 

The rainfall by months during each of the sixty-three 13-month 
periods, between August 1 of the year preceding growth of annual 
ring and August 31 of the year of growth, was expressed by an orthogo- 
nal polynomial of sixth degree. 

In order to determine whether progressive changes in rainfall 
distribution within the 13-month periods could be detected over the 
course of the 63 years under investigation, each of the coefficients 
was plotted on the year of its occurrence, and to each set of 63 values 
an orthogonal polynomial of degree 6 was, in turn, fitted. It was 
found that the a, coefficient, which represents the average monthly 
rainfall during each 13-month period, was the only one which varied 
systematically over the 63 years. 

As a means of eliminating possible progressive changes in annual 
ring width not associated with progressive changes in average monthly 
rainfall, orthogonal polynomials of several degrees were also fitted 
to the ring-width data over the 63 years, and the correlation between 
residuals from trend of these and residuals from trend of various 
degrees of the a, rainfall coefficients were calculated. While the 
highest correlation coefficient existed between residuals from trend of 
degree 3 for annual ring width and from trend of degree 2 for the a, 
coefficient, it was very little higher than the correlation coefficient 
between annual ring width residual from trend of degree 3 and the a, 
coefficient residual from trend of degree zero. As a consequence it 
seemed that the most appropriate expression of the effect of rainfall 
distribution on annual ring width should be obtained from the regres- 
sion of annual ring width residual from trend of degree 3 on the rainfall 
coefficients from which secular trends, if present, were not eliminated. 
This was accordingly done. For comparison, the regression equation 
of annual ring width, from which trend was not eliminated, on the 
rainfall coefficients was also calculated. 

By means of tables of orthogonal functions these regression equa- 
tions are made to express the average effect of an added centimeter 
of rainfall per month during the 13-month period on width of annual 
ring. 

While the regressions based upon the 13-month period were signifi- 
cant, they did not seem to portray the influence of rainfall with suffi- 
cient completeness. Accordingly, the period length was extended to 
15 months, thereby including the rainfall of June and July of the 
year preceding growth. A new set of rainfall coefficients of orthogonal 
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polynomials was calculated, and the width of annual ring expressed 
in terms of these. The results, showing satisfactory completeness, 
were not assuredly significant. Several reasons are suggested in 
explanation. 
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THE INFLUENCE OF RAINFALL UPON TUFT AREA AND 
HEIGHT GROWTH OF THREE SEMIDESERT RANGE 
GRASSES IN SOUTHERN ARIZONA' 


By Pavut B. Lister, formerly assistant forest ecologist, Southwestern Forest and 
Range Experiment Station,? and Francis X. SCHUMACHER, senior silviculturist, 
Forest Service, United States Department of Agriculture 


INTRODUCTION 


It is common observation that the forage stand of the semidesert 
grassland ranges of the Southwest varies greatly from year to year. 
While this variation is undoubtedly due to a number of factors, it is 
generally assumed to be governed predominantly by moisture made 
available through rain—an assumption that is scarcely to be ques- 
tioned. Not all of the rainfall, however, can contribute effectively 
to growth of the forage species. An appreciable proportion of the 
water from showers may run off quickly in certain seasons of the 
year when rainfall is intense; or in other seasons, may evaporate 
directly from the soil in the face of dry winds or high temperatures; 
or the rain may be insufficient to saturate topsoils. 

On the other hand, species or genera in the same plant association 
may not respond equally to the rainfall as the seasons progress; for 
it seems not unlikely that should one species or genus react to the 
same rainfall and in the same way as another species or genus with 
which it is associated, these two must be actively competing for 
moisture, and one of them—that one which is dominant in other 
physiological processes—may eventually suppress and finally supplant 
the other. But should one be particularly sensitive to, say, late 
winter and early spring rainfall while the other benefits largely from 
fall precipitation, they may, if other physiological processes are not 
decisive, grow together in apparent equilibrium. In short, differing 
degrees of response of the associated forage species to a unit of rain 
as the season advances may be a necessary condition of a climax, or 
even of a stable plant association. 

It is the purpose of this paper to report the method and results of 
an investigation into the responses of three important forage grasses 
to variations in the amount and distribution of seasonal moisture 
supply on areas from which grazing animals had been excluded. 
These areas are on the Santa Rita experimental range of the South- 
western Forest and Range Experiment Station in southern Arizona. 


THE DATA 


The climate of the Santa Rita experimental range is characterized 
by long, hot summers with wide variations in temperature between 
day and night; short, mild winters; and low annual precipitation. 

he lowest annual precipitation recorded at the range headquarters, 
4,300 feet elevation, was 12.5 inches in 1903; the highest, 38.2 inches 
1 Received for publication Aug. 26, 1956; issued March 1937. 
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in 1931; and the average from 1901 to 1931, 19.8 inches. Rainfall 
comes in two distinct periods, summer and winter. The summer 
period includes July to September, during which about 56 percent of 
the annual precipitation (July 1 to June 30) occurs; in the winter 
period from November to March about 36 percent occurs. The 
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FicureE 1.—Portion of the Santa Rita experimental range showing location of study plots and adjacent 
precipitation stations. 


remaining 8 percent falls during the more or less dry periods of April, 
May, June, and October. 

Before 1903 the area had been public domain open to uncontrolled 
grazing. In the early summer of 1903, however, the experimental 
range was fenced and ther eafter largely pr otected from grazing 
animals until 1915, when it was again stocked with cattle. 
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Early work was devoted to observation of the natural restoration 
of the forage. In 1915, however, tests of the effect of various degrees 
of controlled grazing were begun. In order to provide check areas 
protected from grazing animals, five enclosures (nos. 3, 4, 5, 6, and 
7 on fig. 1) were established in 1917 and two additional enclosures 
(nos. 17 and 27) in 1925. 

Three precipitation stations, designated by the names of Box, 
Eriopoda, and Huerfano (fig. 1), were located in the vicinity of 
enclosures 3, 6, and 7, respectively ; near enclosures 4 and 5 a precipita- 
tion station designated Rodent was established. Enclosures 17 and 
27 were in the neighborhood of the already established Road precipita- 
tion station. Monthly precipitation was recorded at the stations, 
beginning in the summer of 1922. 

Within the enclosures, on one or more sample quadrats each 1 
square meter, the tuft areas of perennial grasses, 1 inch above ground, 
were mapped with a pantograph (3)* annually at the close of the 
growing season, and measurements were taken of the height of flower 
stalks of the perennial grasses within and around the quadrats. The 
quadrats within enclosures 4 and 5 were not, however, charted afte1 
1929, and occasionally a charting was omitted at other enclosures. 

The effect of rainfall distribution as reflected in yearly changes in 
tuft area and height growth was observed for the three principal 
grasses of the range, viz, three-awn grass (Aristida spp.),* black grama 
(Bouteloua eriopoda), and Rothrock grama (B. rothrockii). From the 
quadrat charts, the difference between tuft areas of the grasses in 
successive years (hereafter called density change) was computed in 
units of square centimeters for each quadrat. In a number of cases 
zero change in density occurred as a result of zero tuft area in each 
of 2 successive years. These were discarded from the study because 
they imposed undue limitation upon a freedom of variation of density 
change which would not have occurred if the quadrats had been large 
enough to preclude temporary extinction of the grasses thereon. 

The density change and height-growth data are presented in table 1. 


TABLE 1.—Measurements of density change and height growth of flower stalks of 
three grasses according to enclosure and year of growth, Santa Rita experimental 
range 

DENSITY CHANGE! 
Grass type and 


. 1923 1924 1925 1926 1927 1928 1929 1930 1931 1932 
enclosure 2 


Square | Square | Square | Square | Squore | Square | Square | Square | Square | Square 


Three-awn: cm em cm cm cm cm cm cm cm cm 
Enclosure 3 (1) — 239 39 —144 ; 2 2 
Enclosure 4 (1) —4Y s —47 —1 
Enclosure 5 (1) —3 —19 —21 
Enclosure 7 (1) —66 2 —19 —5 25 47 68 
Enclosure 17 (1) . ; 44 OF —26 116 65 27 68 
Enclosure 17 (2) ‘ 75 132 —151 —54 1! 16 —4 
Enc.osure 17 (3) 24 20 —18 1] ) 13 —10 
Enclosure 17 (4) . 31 46 —7 51 32 81 26 
Enclosure 27 (2) ‘ 6 6 —2 0 29 —28 —2 

Black grama 
Enclosure 3 (1) —Y 3 —3 P * 3 
Enclosure 4 (1 3A DF howe ae —105 —l4 
Enclosure 5 (1) 36 10 |_. —75 —22 
Enclosure 6. _ _- —18 13 36 —4s8 —8 101 147 —7 





! Averages for all quadrats and years are: Three-awn 6.5; black grama 28.3; Rockrock grama —4.5 
* Quadrat numbers given in parentheses. 
Reference is made by number (italic) to Literature Cited, p. 121. 
4 Mainly Aristida divaricata, but including small amounts, never exceeding 10 percent, of A. orcuttiana, A. 
ternipes, and A. californica 
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TABLE 1.—Measurements of density change and height growth of flower stalks of 
three grasses according to enclosure and year of growth, Santa Rita experimental 
range—Continued 

DENSITY CHANGE—Continued 


Grass type and 


enclocure 1923 1924 1925 1926 1927 1928 1929 1930 1931 1932 
Square | Square | Square | Square | Square | Square | Square | Square | Square | Square 
Black grama—C on. cm cm cm cm cm cm cm cm cm cm 
Enclosure 17 (1) —6 66 —2l 56 8 66 —66 
Enclosure 17 (2) 33 107 —1 ill 38 —2 -Y 
Enclosure 17 (3) —3 73 —33 YS —30 48 54 
Enclosure 27 (1 i i A 31 235 —57 116 —3 Os 122 
Enclosure 27 (2) ‘ : 35 105 —44 107 108 49 69 
Rothrock giama | 
Enclosure 3 (1) 0 —l 20 —171 99 306 49 —60 
Enclosure 4 (1) —64 85 —66 —I " 
Enclosure 5 (1) —3l —30 —6 —132 5 a 
Enclosure 6 . 29 —55 —42 : -| —140 64 223 —70 —114 
Enclosure 7 (1) 109 — 256 95 x Ys —62 85 —170 38 
Enclosure 17 (1). Seetdtie 100 83 | —229 —4 —5 —6 2 
Enclosure 17 (2) = . 22 —-23 | —87 -4 2 —1 —7 
Enclosure 17 (3) 7 pa —Il1 78 —k9 17 68 —94 y 
Enclosure 17 (4) 119 194 —381 37 85 —3 —62 
Enclosure 27 (1) —24 138 — 184 243 15 — 204 23 
Enclosure 27 (2) 38 176 — 294 251 —48 —44 24 
HEIGHT GROWTH 3 
Three-awn Cm cm Cm cm cm cm cm cm cm Cm 
Enclosure 3 (1) 60 50 50 75 ads 45 —_ 8. 72 
Enclosure 4 (1)_- 70 35 P 70 41 70 
Enclosure 5 (1) 65 45 | 36 45 e = 
Enclosure 6 70 50 50 68 | 45 62 85 R5 
Enclosure 7 (1 65 50 76 |}... 75 50 50 65 &5 &5 
Enclosure 17 (1). 45 60 40 55 45 70 70 
Enclosure 17 (2) 45 60 40 55 45 70 70 
Enclosure 17 (3) 45 60 50 | 60 45 R5 76 
Enclosure 17 (4) 45 60 50 60 45 &5 76 
Enclosure 27 (1) ; 35 26 60 46 80 75 
Enclosure 27 (2) , S 35 rae 26 60 46 sO 75 
Black grama 
Enclosure 3 (1)... 55 25 45 60 42 40 52 60 60 
Enclosure 4 (1). 65 35 50 46 50 
Enclosure 5 (1) 60 20 42 | 41 42 
Enclosure 6. . 65 30 55 42 | 40 41 50 50 
Enclosure 7 (1) 60 bee wn 60 48 50 54 53 52 
Enclosure 17 (1) 45 55 40 50 50 63 60 
Enclosure 17 (2) ee 55 40 50 50 63 60 
Enclosure 17 (3) 45 ‘ 45 35 £0 50 65 60 
Enclosure 17 (4) 45 45 | 35 50 50 65 60 
Enclosure 27 (1 , 50 45 50 35 45 48 5S 56 
Enclosure 27 (2) 50 45 50 35 45 48 58 56 
Rothrock grama 
Enclosure 3 (1) 50 35 40 50 36 35 45 60 50 
Enclosure 4 (1) 60 30 40 39 40 
Enclosure 5 (1 55 30 40 21 40 F 
Enclosure 6 55 35 50 38 35 40) 53 50 
Enclosure 7 (1) 55 60 57 42 40 50 58 52 
Enclosure 17 (1) 45 45 35 35 45 58 52 
Enclosure 17 (2) 45 45 35 35 45 58 52 
Enclosure 17 (3 45 7 45 35 25 40 60 50 
Fnelosure 17 (4) 45 45 35 25 40 60 50 
Enclosure 27 (1) 40) 40 50 36 45 40 58 54 
Enclosure 27 (2) 40 40) 50 36 45 40 58 54 


+ Averages for all quadrats and years are: Three-awn 58.7; black grama 49.2; Rockrock grama 44.6. 
METHOD OF ANALYSIS 


The method of analysis was first suggested by the researches of 
Fisher (2) concerning the influence of precipitation on the yields of 
wheat in England. Tippett (8) used the same method in an investi- 
gation of the effect of sunshine on wheat yields. It was used, also, 
in a recent analysis by Schumacher and Meyer (6) of data gathered 
by Dr. Hermann Knuchel in Switzerland relating to the influence 
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of precipitation on the growth of Abies alba. In brief outline, and 
taking density change of grass for illustration, the method is as 
follows: 

The average effect of an additional inch of rainfall above the mean 
during a time interval (¢) within a given period of time on density 
change (d) of a grass between one fall and the next may be expressed 
by the equation 


d= a+ aX, (t) + apX2(t) + -- _-+a,X;,(t) (1) 


in which X,(t), X,(t),_.....X;,(¢) are orthogonal functions of the time 
interval (t) during the period, and ap, a, a2,-.-.--a,, are constants 
to be determined. 


If the rainfall (7) within the same time interval (¢) of the same period 
may be expressed 


r=d)+a,X,(t)+a.X2(t)+__---- +a,X;,(t) (2) 


then the whole effect of the rainfall distribution during the period on 
change of grass density (D) from one fall to the next may be written 


T T 
t=— t= 5) 
D—Mp= > dr — > (a+ aX, (t) + Q.X9(t) TT o2esam 
T T 
t=-; t=—— 
2 oy, X ,.(t)) (ap + aX, (t) + a2.X0(t) + _- . +a,X;,(t)) (3) 


where /, is the mean value of the observed annual changes in grass 

density. The summation is carried over all time intervals of the 

period, the first and last of which are denoted by —3 

tively, 7 being the number of time intervals within the period. On 

account of certain properties of the orthogonal functions X;,(¢)[? 
1 


7 
and +,;, respec- 


=, 1,2 kj, such as [X,=0; T ~X"*,;=1; and >X,X,;=0 where 
i ~ j, [j=0, 1, 2,__-___k], equation (3) becomes 


D— Mp=T (apo + a 0,+ and.+ ____-- + a,0;) (4) 

The object of the analysis is the pictorial representation of equa- 
tion (1). The coefficients a), a,, d:,-.....a, of equation (2) are 
first determined by the method of least squares from the rainfall 
record of each period at the precipitation stations. As many obser- 
vation equations of the form of equation (4) are then set up as there 
are values of annual density change to be paired with the rainfall 
coefficients. Normal equations are derived from the observation 
equations according to the method of least squares and their solution 
affords estimates of the unknowns ap, a, a,---.--a,. Finally, 
when these are put in equation (1) the average effect appears of an 
additional inch of rain in each time interval on the annual change in 
grass density. 

Period lengths of 15 months were chosen, beginning with June 1 
of the year preceding growth and ending with August 31 of the 
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current year of growth. While September rainfall of the current 
year may also play a part in grass growth, it happened that during 
2 years the quadrats were charted in September before the last of the 
summer rains. The effect of September rain on the current year’s 
growth might not, as a consequence, be adequately evaluated, and 
it was therefore dropped from consideration. 

The extension of the initial limit of the period beyond September 
was suggested by the work of Nelson (5), who pointed out that, in 
black grama range in southern New Mexico, the rainfall of the pre- 
ceding summer effects significant changes in plant area from one fall 
to the next. He concluded that the value of the rainfall during the 
preceding summer was carried over into the current year’s growth 
through the vigor imparted to the growing plant. Accordingly, to 
include the relatively high summer precipitation following compara- 
tively dry spring weather, June 1 was chosen as the initial limit of the 
period. With the rainfall record by months, thé time interval ¢ is 1 
month, and 7’, the number of intervals in the period, is 15. 


THE RAINFALL DISTRIBUTION CONSTANTS 


As stated above, the first part of the work is the determination of 
the coefficients dp», a, ds, _a, of the orthogonal polynomial 
equation (2) for each period for which data are at hand. An imme- 
diate consideration is the degree of polynomial needed in order to 
describe adequately the rainfall distribution in each of the periods. 
Fisher (2), using the rainfall record at Rothamsted, England, cal- 
culated coefficients up to degree 5, but his data were according to 
intervals of 6 days in a period length of 1 year. Schumacher and 
Meyer (6) used sixth-degree polynomials to describe the precipitation 
distribution at Aarau, Switzerland, for time intervals of 1 month in a 
period of 15 months. Precipitation in the semidesert ranges of the 
Southwest is, however, relatively much more variable than in locali- 
ties of greater abundance of rain. In view of this variability and the 
relatively long time interval employed, polynomials of degree 8 were 
decided upon. 

The computational work, though tedious, is not complex; for if there 
are equal intervals of the independent variable, and if to every value 
of the independent variable there corresponds only one value of the 
dependent variable—as rainfall per month—all the computations 
necessary to calculate the curve representing the trend of (r) with 
respect to (¢) for each period may be decidedly simplified through the 
use of orthogonal functions, as these have been tabulated by different 
authors. A very convenient set of tables is given by Lorenz (4), 
together with a full discussion of the theory of orthogonal functions 
and their application to specific problems. 

Following the computational scheme and checking method as recom- 
mended by Lorenz—known as the scheme of Chetverikoff—the rainfall 
coefficients dp, a), @2,....--ds, of equation (2) were calculated from the 
rainfall record of each period at each station. These are listed in table 
2 for units of 0.01 inch of rain per month. The rainfall trend for the 
average of the periods and stations as derived from the arithmetic 
average of each coefficient, is compared with the recorded average 
rainfall in figure 2. 
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FIGURE 2.—Comparison of curve of rainfall trend (derived from arithmetic average of each coefficient) with 
recorded average monthly rainfall. 


TABLE 2.—Coefficients of orthogonal polynomials which express precipitation trends 
(in hundredths of inches per month) from June of a preceding year to August 
inclusive of a current year, according to precipitation stations; Santa Rita experi- 
mental range 


Rainfall coefficients ! 








™ ; Cur- 
Precipitation sta- rent 
tion 
vear 
¢ a a ( as Gs (6 a as 
124 25 97 67 35 16 70 —53 —40 
153 ~39 35 35 34 47 89 31 —82 
sy 24 4 49 —2 42 36 25 —21 
126 16 48 &5 29 52 —1| 3 —5 
Box... 122 —14 42 80 —12 3 31 —25 45 
116 18 51 62 7 8 -31 —17 3 
155 —4 32 37 —56 48 31 —50 41 
207 107 156 126 112 75 5 4 —22 
248 —52 116 109 —54 69 —73 —99 76 
151 35 121 80 —36 | 48 —85 —59 —23 
| | 152 —64 41 35 — 50 75 —S6 35 —68 
Rodent . : 121 —15 24 71 —14 36 —8 —25 4 
102 -9 30 56 3 29 34 —20 19 
95 —2 28 41 —11 39 ~30 —21 1 
132 5 | 98 67 —41 53 —81 —25 Zi 
148 —2 37 34 —44 23 | —100 —1 —79 
100 4 77 | 43 —2 39 —56 51 —33 
128 14 50 90 18 | 50 4 —5 —15 
Eriopoda 119 —-2B 35 69 —32 | 40 —42 — 50 36 
| 115 15 54 65 12 56 29 —8 —15 
171 26 | 50 | 46 —63 | 34 — 59 —71 28 
204 63 | 164 | 107 89 | 92 —27 57 —5 
207 —45 100 | 91 | —63 52 —93 —94 4 
142 11 | 114 | 81| —24 60 —69 —20 —32 
161 —25 37 | 44) —42| 35 —95 0 —70 
100 16 SY 414 5 53 —53 54 —23 
158 —10 48 89 | —1 |} 60 —11 1 —24 
Huerfano 142 —37 36 90 —40 | 46 —19 —58 4s 
125 14 5l 67 —4 46 —37 —36 —13 
170 18 53 47 —45 58 —51 —37 29 
217 73 157 116 94 | 85 —10 sy 39 
230 —60 95 90 —65 | 40 —88 | —115 49 
87 24 58 37 | —12 14 —34 14 —5 
148 —20 37 30; —3| 26 =] —32 —23 
149 —13 40 82 —16 | 50 —10 —29 9 
Road 122 —8 49 75 —1 | 45 —33 — 2% 25 
, 129 9 | 47 52 —31 | 36 —49 —39 5 
181 —2) 48 17 —49 | 42 —27 —40 33 
210 89 173 W9 103 | 58 —5 13 14 
230 —64 106 102; —54 82 —78 —63 52 





1 Average ao, 149.425, a;, 1.925, a2, 70.200, a3, 68.025, a4, — 11.600, a5, 49.050, as, — 43.725, az, — 18.675, as, — 1.000 
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REGRESSION OF ANNUAL DENSITY CHANGE AND HEIGHT GROWTH 
ON RAINFALL IN DIFFERENT SEASONS 


The unknowns a, a, a, --.. a of equation (4) were next calculated 
by the method of least squares (1, ch. XJV) from the data of tables 1 
and 2, after setting up observation equations of the dependent variable 
(density change or height growth of one of the grasses, as the case might 
be) in terms of the rainfall coefficients ay, a), d2, _... as for correspond- 
ing station and year. As an illustration, the observation equation for 
change of density of three-awn grass between the falls of 1922 and 1923 
at Box station is 


—245.5=124Ta+25Ta,+ _- —40Ta,+C 


where C is a constant, also to be determined. Each set of observation 
equations leads to the normal equations which, when solved simulta- 
neously, render 15 (since 7 =15) times the coefficients Qip, Ot, Qa, ---- As 
of equation (4) and the constant C. Their numerical values are listed 
in table 3 according to the density change or height-growth equation 
to which they belong for each of the three grasses. 


TABLE 3.—Regression coefficients and constant terms of equations which express 
calculated departures of density change or height growth of flower stalks of 3 
range grasses from their respective averages, in terms of the rainfall coefficients of 
the 15-month periods 


DENSITY CHANGE 


Regression coefficients 





Species - ——- ~ 7 " ( onstant 
term 

5a, L5a l5a l5a l5a l5a ber. l5ax i See 

2 ‘ B) 4 7 8 i] 
Three-awn 1. 1349] 0. 4985] —1. 5968/—0. 2768; 0.7292) 1.8153) 0.0467| —0. 7048) —9. 4643) — 139. 7243 
Black grama 1.4772) 2.0174|—1. 3738) 2.6970|—2. 5859) — 3¢ 1.5715) 1.5045) —. 6613) —217. 3320 
Rothrock grama__- 4456) 1. 5617 . 6979 . 7796| —4. 2941) 1 4.0215 . 1654) —1. 9016 —125. 4921 

HEIGHT GROWTH 

Three-awn 0. 1284) 0.2078) 0.0404) 0. 3803) —0. 2269) —0. 0378) —('. 0226) 0.1102) 0.0028) —47. 7886 
Black grama .0156; . 0964 . 2639 . 1367) —. 3251 . 0124 . 3295) —. 0258) —. 0854) —21. 6410 
Rothrock grama 0147| =. 0885 2489 1490) —. 2654) —. 0533 . 2442 0672; —. 0639; —18. 7090 


Tests of significance are based upon analyses of variance 
presented in table 4. In general 


~_~ 


7) as 


~(y—M,)?=z(y— Y)?+2(Y—M,) 


in which y is an observed value of the dependent variable, the arith- 
metic average of which is M,, and Y is the estimate of y by the 
method of least squares. The left-hand member of the above identity 
is the total sum of squares of density change or height growth before 
eliminating that part due to their regression on rainfall distribution 
during the 15-month periods. The sum of squares due to regression 
is calculated through the identity 


=(Y¥—M,)?= TaoX (agy) + Tayd (ayy) + Tad (azy) + ~~~. + Tagd (agy) 


for which the sum of products of each term of the right-hand member 
was computed previously as steps leading to the normal equations. 
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The remaining sum of squares, 2(y—Y)*, is due to deviation from 
regression, that is, it is that part of the total sum of squares that has 


been found to be independent of the rainfall distribution. 


TABLE 4.- 


Range gruss and source 


rhree-awn 
Regression __._- 
Deviation from regression 


Total_- 
Black grama: 
Regression ---- 
Deviation from regression 


Total... 


Rothrock grama: 
Regression 


DENSITY CHANGES 


2 of variation 


Deviation from regression _ -- 


Total___- 


Three-awn 
Regression 
Deviation from regression 
Total 


Black grama: 
Regression 


Deviation from regression -_ - 


Total. 
Rothrock grama: 
Regression... - 
Deviation from regression 


Total 


HEIGHT 


Sum of | Degrees of 


squares | freedom 

71, 241 y 

122, 339 44 

193, 580 53 

103, 318 | 9 

116, 852 | 45 

220, 170 | 54 

397, 988 i) 

676, 255 64 

-|1, 074, 243 73 
GROWTH 

11, 664 9 

6, 077 63 

17, 741 72 

6, 368 9 

2, 420 69 

8, 788 78 

4, 878 9 

1, 683 69 

6, 561 78 


Mean 
squares 


7, 916 


2, 780 


11, 480 


2, 597 


44, 221 
10, 566 


707.6 
35. 1 


542.0 


-Analyses of variance of density changes and height growth 


F ob- F where 


served | P=0.01 


| 


2. 85 2. 85 
| 

4.42 | 2. 84 
| 

4.19 2. 70 

13. 43 2.71 

20. 16 2. 68 

22. 21 2. 68 


Sums of squares due to these three sources of variation are listed 
in the second column of table 4. 
of freedom—one less than the numbers of observations for total, one 


for each independent variable (a, a,, a2, —- 
and the remaining number for departures from regression. 


fourth column are the mean squares. 

If density change or height growth of the grasses were independent 
of the rainfall distribution, the mean square due to regression would 
be approximately equal to the mean square of the departures from 
regression, and the ratio of the former to the latter—designated F- 


would therefore be 1, approximately. 


In the third column are the degrees 


_. ds) for the regression, 


In the 


Under such conditions, values 


taken from the distribution of Fin random samples (7, pp. 88-91) as 
given in the last column would be exceeded but once in 100 sets. By 
comparison, the observed values of F in the next to last column are 
All the correlations of grass growth with the 
rainfall distribution are consequently adjudged highly significant. 
The issue, however, is not the obvious truth that rainfall influences 


as great or greater. 


grass growth; but since the data at hand are sufficient to corroborate 
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what is already known, the issue is, rather, to indicate the extent of 
rainfall influence at different seasons of the period. 

Putting, therefore, the set of coefficients ap, a, a2, ---- as, which 
express annual change in density, or height growth of one of the 
grasses concerned, into equation (1), where 


d= a+ aX, (t) + a2X2(t)+----+asXs(t), 


the functions X,(t), X2(t), X,(t) being taken from Lorenz’s 
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FIGURE 3.—Computed effect of an additional inch of rain per month on annual change in density of three 
I I 
grasses. 
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FIGURE 4.—Computed effect on yearly growth of flower stalks of an inch more of rain per month. 


tables (4),°> and plotting on ft, we have the graphic expression of the 
effect of an additional one-hundredth inch of rain per month, accord- 
ing to the months of the period, on annual change in density or height 





5 Lorenz’s tables do not give the functions X;(¢ for i greater than 6. Consequently X7;(¢) and Xs(f) were 
calculated for each of the 15 values of ¢ through the use of equation (11) of Lorenz's paper (4, p. 3/). 
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growth. The results (multiplied by 100) are presented graphically in 
figure 3 to show the average effect of 1 additional inch of rain per 
month on annual density change in square centimeters per meter 
quadrat; and in figure 4, on height growth of flower stalks in centi- 
meters per year. The ordinates of these figures represent deviations 
from the mean values of density change or height growth listed at 
the bottom of table 1 


DISCUSSION 


The curves of figures 3 and 4 represent the average benefit to the 
forage stand, in tuft area and height growth, of 1 inch of rain per 
month during the 15-month period in excess of each monthly rainfall 
average as given by the distribution curve of figure 2. Change in 
sign of the ordinate scales would, of course, represent the correspond- 
ing average damage chargeable to 1 inch of rain per month in defect 
of its monthly average.® Likewise the best estimate of density change 
or - ight growth as a result of no change from the monthly averages— 


i. , the average rainfall distribution ‘itself —would ie zero deviation 
rio the averages given at the bottom of table 1—1i. e., the general 


averages themselves for all stations and years. 

It is the experience of range workers that the three grasses under 
discussion associate in apparent tolerance of one another. The trends 
of figure 3 suggest that this may be due to differences in response of 
the grasses to rainfall distribution. The first general impression 
conveyed by the figure is that rainfall does not affect the density 
changes of the three grasses alike, except at seasons—late July of the 
preceding and current years—when the effect of rainfall variation is 
rather negligible. 

Closer inspection of figure 3 reveals that the curves of the grama 
grasses show approximately the same trends, the differences between 
them being that black grama lags by about a month behind Rothrock 
grama and is, moreover, less sensitive to rainfall variation of the 
winter and spring. 

The average tuft area per meter quadrat during the years of obser- 
vation was 168 cm? of black grama, 162 cm? of Rothrock grama, and 
76 cm? of three-awn. For reasons that are not clear, additional 
rainfall during June of the preceding year results in loss to the current 
year’s density growth of black grama at an average rate of 76 cm? 
per inch (45 percent of average tuft area); of Rothrock grama at an 
average rate of 53 cm? per inch (33 percent). Although variation of 
July precipitation of the preceding year has little effect on density 
change of either of the gramas, additional rain during the late summer 
benefits both. Maximum benefits, coming from the late September 
excess to Rothrock grama and from the late October excess to black 
grama, reach 56 cm? per inch or somewhat over 30 percent. 

Additional rainfall during the winter season, on the other hand, has 
a deleterious influence on the density growth of the gramas, particu- 
larly on that of Rothrock grama, which is subject to an average loss 
of 80 em? (49 percent) per inch of additional January rain. The 
greatest loss to black grama density is 26 em? (15 percent) per inch 
of additional February rain. 

Extra spring precipitation is beneficial to the gramas, its maximum 
effect reaching 70 cm? (42 percent) per inch of April rain to Rothrock 


6 In this case it is recognized, of course, that as much as an inch in defect of most of the monthly averages 
from October to June is impossible 
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grama and 44 em? (27 percent) per inch of May rain to black grama. 
Thereafter the benefit is lessened, and in middle summer it is negligible. 
August excess, however, again benefits the density growth of the 
gramas to an extent of about 30 cm? (18 percent) per additional inch. 

The reaction of three-awn tuft area to precipitation is quite distinct 
frcm that of the grama grasses. Except during June and early July, 
additional precipitation during the preceding summer and fall has 
little effect on its density. An average density loss of 20 em? (26 
percent) is associated with an additional inch of June precipitation in 
both the preceding and current years. 

During the winter season, additional rain becomes progressively 
more beneficial to the tuft area of three-awn, reaching its maximum 
in March with 45 em? (59 percent) per inch. Additional April pre- 
cipitation is less beneficial, that of May ineffective, and that of June, 
again, depressive. Thereafter tuft area of three-awn is but negligibly 
affected by variation in rainfall. 

Turning now to figure 4, it is at once apparent that additional rain 
coming during the preceding fall and current spring and summer bene- 
fits the height growth of the three grasses during the summer of the 
“current year”, but coming during the winter it mostly detracts from 
it. The practically identical effect of the rainfall distribution on the 
height growth of the two gramas, in contrast te its singular influence 
on three-awn, is also striking. 

During the years of observation the average height of the flower 
stalks of black grama was 49 cm; of Rothrock grama, 45 cm; and of 
three-awn, 59 em. 

Until late summer of the preceding year, additional rainfall has 
little influence on the current year’s height growth of the grama flower 
stalks. But with the advent of August the efficiency of additional 
rainfall mounts, reaching a maximum benefit in late September of 
5 em per inch for black grama and 4 em per inch for Rothrock grama, 
or 9 to 10 percent of the average height in either case. As fall ad- 
vances, the benefit lessens, dropping through zero in November to a 
loss as a result of excess winter precipitation, which reaches 7.5 cm 
per inch for black grama and 6.4 em per inch for Rothrock grama (14 
to 15 percent) in January. Additional spring rainfall is again bene- 
ficial, particularly in May, when the benefit reaches 3.5 em per inch 
for black grama and 3 cm per inch for Rothrock grama (7 percent). 
Thereafter the benefit is lessened to a minimum of about 15 em per 
inch in July before increasing again with advancing summer. This 
same effect of July rain on density change is apparent in figure 3. A 
not unlikely reason for this seeming loss of rainfall efficiency is that 
with the bursting of flower buds the benefit of additional rain is with- 
drawn from height and tuft area growth and is put into flowers and 
seed. 

Roughly, the rainfall distribution affects height growth and density 
change of the two grama grasses alike, for the curves of figures 3 and 4 
show generally the same trends. In three-awn, on the other hand, 
height of flower stalks and density change are affected quite differently 
by additional rain. For whereas additional rain during the preceding 
summer and fall scarcely affects density change at all, it contributes 
at an average rate, from August to the end of the preceding year, of 


2.5 em (4 percent) per inch per month to height growth during the 
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summer of the current year. During the later winter, when additional 
rain greatly benefits tuft area, it detracts from height growth at a rate 
that reaches 3 cm (5 percent) per inch of additional March rain. By 
June the influence of rain is again reversed, an additional inch effecting 
a loss of 26 percent to density but a benefit of about 3.5 cm (6 percent) 
to height growth. Through July and August this benefit is somewhat 
increased. 
SUMMARY 


The effect of variation in precipitation distribution during a 15- 
month period on the annual density change and height growth of the 
three forage grasses—three-awn grass (Aristida spp.), black grama 
(Bouteloua ervopoda), and Rothrock grama (B. rothrockii)—within the 
growing season of the current year, was determined by statistical 
analysis. The method used consisted in describing the observed 
monthly precipitation during each 15-month period, ending with 
August for each station and year, by orthogonal polynomials of 
degree 8; the coefficients of these were then used as independent varia- 
bles in six multiple-correlation analyses, the dependent variables of 
which were (1) the annual density change and (2) the height growth 
of flower stalks, of each of the three grasses according to station and 
year. All the correlations are highly significant. 

Graphic expressions of the average effect of an added inch of pre- 
cipitation per month indicate, in the case of the two grama grasses, 
that the most beneficial precipitation distribution from normal for 
both density change and height growth consists of relatively dry win- 
ters sandwiched between relatively wet autumns and springs. ‘Three- 
awn grass seems admirably adapted to the variation in precipitation 
distribution characteristic of the area studied in that either positive 
or negative departure from normal seasonal precipitation benefits, in 
general, either density change or height growth, although not both 
at once. 
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INTRODUCTION 


To be able to estimate accurately the number of wireworms in any 
given field is a fundamental necessity if progress is to be made in 
economic investigations of these insects. With valid methods for 
estimating the size of wireworm infestations in the soil the following 
important objectives are attainable: (1) The relative effectiveness of 
various practices for controlling wireworms can be determined by 
statistical treatment; (2) the actual degree of crop damage inflicted 
by different numbers of wireworms can be measured; (3) knowledge 
of the changes in degree of infestation in the field, annually or over 
longer periods, in conjunction with a study of soil factors, will help 
toward the discovery and measurement of the factors affecting wire- 
worm incidence. 

None of the papers on quantitative wireworm investigations thus 
far published refer to the degree of accuracy obtained by the sampling 
method used. A review of the literature indicates that sampling 
methods for censuses of wireworm populations should be studied to 
determine which method gives the most accurate results. 

In studying wireworm populations several methods of sampling 
were tried. At first the samples taken * were large, the units being 
| square yard or more in surface area. These large units were too few 
in number, and the labor involved in their examination proved exces- 
sive. Efforts were therefore directed toward developing a time- and 
labor-saving system that would give fairly accurate results. Smaller 
sampling units, of 1 square foot and less, were studied and will be 
reported upon in this paper. 


REVIEW OF LITERATURE 


Several British investigators refer to their soil sampling of unit 
blocks 9 inches square on the surface and 9 to 12 inches deep. Thus, 
Roebuck (9)* dug out portions of soil in different parts of a field, 
carefully crumbled and sifted it, and counted the wireworms. The 
size of his sampling unit was at first a 1-foot cube, but later a block 
9 inches square on the surface and 1 foot in depth was used. The 
number of units in a field ranged from 5 to 12, and in subsequent counts 
on the same field the same number of units were taken. The average 
number of wireworms found in the unit samples from a field was then 
multiplied b by a conversion factor to obtain an estimate of the number 





1 Received for © public ation June 18, 1936; issued March 1937. 
? The work reported in this paper was done under the direction of M. C. Lane, entomologist, in charge of 
Pa acific Northwest a investigations, Walla Walla, Wash. M.C. Lane, F.H. Shirck, R. S. Lehman, 
and K. E. Gibson assisted in the accumulation of the data 
’ Wireworm-population studies were begun in 1926 by M. C +. and F. H. Shirck at Toppenish, Wash. 
4 Reference is made by number (italic) to Literature Cited, p. 1: 
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per acre. The sampling units of both Morris (7) and Edwards (1) 
were 9 inches square on the surface and 9 inches deep. In his faunal 
studies of arable lands at Rothamsted, Morris examined as many as 
23 sampling units from each half-acre plot used, and states that they 
were never taken nearer together than about a yard. The unit samples 
were carried to the laboratory and, after being dried on brown paper, 
the insects were counted in a small quantity of soil at a time. Ed- 
wards examined an average of four units per plot in his insect surveys. 
In sampling for wireworms King and Glen (4) and Wardle (11) both 
used units 1 square foot in area, King and Glen taking 10 units per 
plot. 

In an earlier publication King (3) refers to the use of a smaller unit, 
0.1 square foot, but he does not say how many of these units were 
taken. For preliminary investigations Hawkins (2) used units 1 
square foot in surface area, extending to a depth of 6 to 8 inches, but 
he does not state how many units were taken. Lacroix (5) determined 
the size of larval populations by excavating two areas, each 2 feet 
square, to a depth of 2 feet. 


METHODS AND PROCEDURE 


APPARATUS AND EQUIPMENT 


The difficulty of sampling is increased because, since the wireworms 
live below the surface of the ground, there is no indication of their 
whereabouts. As there is no natural unit of measurement, arbitrary 
units had to be adopted. In this study the standard sampling unit 
comprises the soil under an area 1 square foot in size, the depth ranging 
from 1 to 2 feet according to the season. The unit will therefore be 
designated by area, it being understood that the depth chosen was 
considered sufficient for obtaining all wireworms. The 1-square-foot 
unit is large enough to make it certain that wireworms will be obtained 
in the majority of cases, and it is as small as can be used where depth 
records are desired for studies of biology and control. 

For taking the samples a metal form of 10-gage black iron, 1 foot 
square and 8 inches high with a 1-inch flange at the top, was devised 
(fig. 1, A). This form is pushed into the ground and the soil removed 
in 3-inch layers with a square-edged shovel. Each layer approximately 
fills a 10-quart galvanized-iron pail, which makes an easy measure. 
The number of pails of soil per unit area depends on the depth of the 
sampling. 

Two smaller sampling units were also investigated. One of these 
has a surface area of one-fourth square foot. A metal form of 8-gage 
black iron, 6 inches square at the bottom, 6% inches square at the top, 
and 12 inches high, was made with a 1-inch flange of 1-gage iron welded 
to the top edge (fig.1,B). Atamper was made of a steel plate 84 inches 
square and five-eighths of an inch thick, with a chromium-steel auto- 
mobile drive shaft (1% inches by 4% feet) riveted loosely at one end 
through the center of the plate (fig. 1, D). The form is driven into the 
soil its entire depth with the tamper, and then the 4-cubic-foot section 
of soil is removed intact by using the steel plate of the tamper under 
the flange asa pry. The soil is then shaken into a pail for sifting. The 
form, being slightly larger at the top, allows easy removal of the soil 
by shaking. Soil below a depth of 12 inches is removed with a soil 
auger (fig. 1, F). 
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The smallest sampling unit studied is one-sixteenth of a square foot. 
The same type of square metal form was used as with the larger units, 
the lower opening being 3 inches square and the top opening 3% inches 
square to allow easy removal of the soil (fig. 1, C). This size is the 
most convenient of the three to use, as the metal form is more easily 
driven into the soil and lighter to handle. The unit of soil is small 
enough to be run through a washing device (fig. 2) (10), whereby all 
stages of wireworms, even the eggs, can be almost completely 
recovered. 

Very early in the investigations a mechanical soil sifter was developed 
to reduce labor and increase the number of sampling units possible to 
be taken in a given time. The first machine (6) devised was operated 
by hand. A more mobile power sifter * has finally been developed with 
sifting capacity of approximately 5 cubic feet per hour, although the 
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FIGURE 1.—Equipment and forms used in taking soil samples to determine wireworm populations: A, 
l-square-foot sampler; B, '4-square-foot sampler; C, e-square-foot sampler; D, steel tamper; E, soil 
augers for use in obtaining samples deeper than first surface foot. 


number varies according to the moisture content, texture, and com- 
pactness of the soil. This sifter in its final development is pictured in 
figure 3. For best results a screen of 4 meshes to the inch is used on 
top for breaking up lumps and removing trash. Under that a 10- or 
12-mesh screen can be used to catch adults and large wireworms. 
The mesh of the bottom screen should be small enough to catch the 
smaller wireworms yet large enough to allow most of the soil to pass 
through. The usual size is 16 mesh, but it is safer to employ 18 mesh 
with dry soils. Experience has demonstrated that never more than 10 
percent of the small wireworms (5-10 mm) pass through the*sifter 
screens, the loss being less when the soil is wet. When it is desired to 
obtain census counts of eggs and very small or newly hatched larvae, 
the soil should be washed through screens of 40 to 60 mesh. The most 


5 LANE, M. C., and Suirck, F. H. A MOBILE POWER SOIL sirTeR. U.S. Dept. Agr., Bur. Ent. and 
Plant Quarantine Multigraphed Cire. ET-70, 2 pp., illus. 1936 
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FiGuRE 2.—Soil-washing device used in obtaining wireworm populations in Me-square-foot units of soil. 





Ficure 3.—Power soil sifter used in obtaining populations of wireworm adults and larvae, 
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accurate population counts that are obtained by soil sifting are taken 
during the spring months, before the current year’s brood is hatched. 
Only one pail of soil from 1 square foot of surface to a depth of 3 
inches is placed on the sifter at one time. The numbers of wireworms 
found are recorded on standard forms (figs. 4 and 5) according to the 
depth, their stage of development, and the condition of the soil. 


Form £Q-20 Date 





7008 Place —— ~ ———————— Notes by 





SOIL DATA 
Som, LARVAE REMARK 
Depth | Temp Me 


LARVAE PUPAE ADULTS 





Baurie 
No. 




































































FiGURE 4.—Form used for recording soil-sifting data. 


SOIL MOISTURE DATA 


Som Can Tota, Weicat Tota. Weicat Net Weicat Weicat or MOISTURE 
Derta No. Wer Dey Dry Water Pex Cent 





Saumrie 








15°-18° 
18°-21" 
21°-24" 






































FIGURE 5.—Reverse of form shown in figure 4, used for recording soil-moisture data. 
METHODS OF SELECTING UNIT AREAS 


In selecting unit areas it is not possible, strictly speaking, to obtain 
a theoretically random sample, but the plan followed was to take 
representative unit areas in each field. Sampling was begun at one 
end of a field and continued across to the other end. The practice was 
to scatter the unit areas over the field as evenly as possible, each 
succeeding area being located by a toss of the metal form. This plan 
has proved the most feasible and comes nearest to giving a theoretically 
random sample of the wireworm population in a field. 

Because of the difficulty in obtaining samples, it is of scientific value 
to know whether the soil in the smaller unit areas (one-fourth or 
one-sixteenth square foot) will give as good an estimate of wireworm 
density as the soil in the larger unit of 1 square foot. Again, it is 
129083—-37—_—-4 
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important to know whether a sample of 25- or 50-unit areas will be 
as good as a sample of 100-unit areas. The size of the unit area and 
the number of units necessary to constitute a sample are discussed 
in this paper. 

During the past several years data on the wireworm populations 
of a large number of fields have been obtained for this statistical study, 
several thousand unit areas having been examined and the wireworm 
count recorded for each. 


PRESENTATION OF DATA 


DETERMINATION OF THE TYPE OF WIREWORM DISTRIBUTION 


In the present problem the variation in the mean number of wire- 
worms from a given number of sampling units would be given by the 
standard error of that mean. This statistic is valid provided the 
sampling units differ from a homogeneous population solely by errors 
of random sampling. It is necessary, therefore, that a test of the 
homogeneity of the wireworm distribution be made. 

In table 1 are given the observed and the expected frequencies of a 
Poisson series (Pearson (8)) for unit areas of 1, 4%, and \%¢ square foot. 
The graphs representing these frequency distributions are given in 
figure 6. In applying the x? test to such a series, in no class should 
the number expected be less than 5, since the calculated distribution 
of x? is not closely realized for very small classes. 

TABLE 1.—Comparison of the observed and expected frequency distributions of 


wireworm-population samples obtained from 1-, \4-, and Y4o-square-foot unit areas, 
each sample consisting of counts from 100 units 


Frequency of occurrence in sampling units of 


Wireworms per unit area (number) 1 square foot 14 square foot i446 square foot 


Observed | Expected | Observed | Expected | Observed | Expected 


0 43 41 81 82 
l 116 18 40 37 16 16 
2 13 13 Ss 16 23 22 
3 17 19 3g 56 
4 14 19 
5 14 16 
6 q ll 
7 : 6 7 
8 and over 11 7 
Total 100 100 100 100 100 100 
x? 12. 47 5. 84 51 
Degrees of freedom 6 2 1 
Probability . 05 . 06 . 40 


! | or 0 wireworm per unit area. 
22 or more wireworms per unit area. 
3 or more wireworms per unit area. 


In testing the 1-square-foot unit the frequencies for 0 wireworm 
and 1 wireworm and also those for 8 or more wireworms are pooled. 
The comparison obtained is given in table 1. When eight frequency 
classes are used, x’=12.47. In ascertaining the value of n, the ex- 
pected frequencies are calculated, not only from the total number of 
values observed (100), but also from the observed mean. There 























jan. 15,1987 Practical Field Methods of Sampling Soil for Wireworms 129 


remain, therefore, 6 degrees of freedom, and n=6. The x? table 
shows that for this value P is between 0.10 and 6.05, agreeing closely 
with expectation. 

In like manner the smaller sampling units are tested. According 
to the x? test for each of these surface areas, the actual wireworm 
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populations fall into Poisson series, indicating the uniformity of the 
population and the usefulness of the technique employed. In the 
case of the 1-square-foot sampling units, the distribution is moder- 
ately asymmetrical (fig. 6, A), the class frequency decreasing with 
greater rapidity on one side of the maximum than on the other. 
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Because of the preponderance of the zero-infestation units, the other 
two distributions obtained from the small surface areas represent 
reversed J-shaped curves (fig. 6, B, C). 


COMPARISON OF SIZES OF SAMPLING UNITS 


Tests were made in seven fields to compare the 1-square-foot, the 
¥-square-foot, and the \s5-square-foot unit areas. The population 
statistics obtained by the three methods of sampling in the same field 
are presented in table 2. 


TABLE 2.—Comparison of the mean wireworm populations in 1-square-foot, \4- 
square-foot, and »-square-foot unit areas obtained from the same field 


Coeffi- | Relative 


. Unitsin | Sizeof | Average wire- 
Field no. “ “ eel cient of | standard 
sample | unit area worms vasiatinns| Gaver 
Square 
Number feet Number Percent Percent 
j 25 1 5. 40+0. 79 73 14. 6 
l 25 ly - 96+ .24 125 25. 0 
| 25 Vie 36+ .11 152 30. 6 
| 25 1 5.444 .91 M4 16. 7 
2 25 4 1.564 .27 86 17.3 
| 25 Ye 80+ . 21 131 26.3 
| 25 l 8. 72+1.33 76 15. 3 
3 25 14 2.044 .50 122 24.5 
| 25 Li 64+ .17 132 26. 6 
| 25 1 19. 00+3. 20 84 16.8 
4 25 4 5. 5641. 20 108 21.6 
| 25 Me i.52+ .41 135 27.0 
| 5O 1 5.424 .59 77 10.9 
5 50 4 1.264 .20 112 15.9 
| 50 li6 58+ .12 151 20.7 
| 50 1 13. 55+ 1. 80 96 13.3 
6 iO ly 3.732 . 66 126 17.7 
| 50 Lie 110+ .22 144 20.0 
| 100 1 9. 6441. 06 109 11.0 
7 100 4 2.534 .38 150 15.0 
| 100 iF 83+ .13 159 15.7 


!The standard error and not the probable error is used after the + sign. 


In comparing samples consisting of the same number of units, 
i. e., 25, 50, or 100, it will be seen that the coefficient of variation for 
the 1-square-foot unit is consistently smaller than for the other two 
units. The relative standard error, which may be defined as the 
standard error expressed as a percentage of the mean, likewise is 
consistently smaller. By the consistency of results it is shown that 
the 1-square-foot unit gives a better estimate of the wireworm popula- 
tion than the smaller units. The \-square-foot samples are also 
shown to be more accurate than the \\s-square-foot samples. 


NUMBER OF UNITS NECESSARY TO CONSTITUTE A SAMPLE 


One-square-foot, one-quarter-square-foot, and one-sixteenth-square- 
foot units were taken in bare fields and in crop-covered fields ranging 
in size from a small fraction of an acre to several acres. In table 3 are 
given the mean population levels for 25 fields, obtained by taking 25, 50, 
and 100 units of each size from the same field, together with their stand- 
ard errors and the relative standard errors. For the smaller units the 
population means, naturally, are smaller. These figures can be placed 
on a comparable basis by multiplying by the proper conversion figure, 
4 or 16. An explanation might be made at this point regarding the 
very large number of wireworms obtained in some of the fields by the 
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‘ie-Square-foot units. In these cases samples were taken in fields to 
determine the total population of newly hatched wireworms. The 
statistics in table 3 show that the relative standard error is not only 
smallest but most stable for the 100 units. The 50 units show good 
results also. 


TABLE 3.— Mean wireworm-population levels obtained by talking samples of 25, 50, 
and 100 units of 1 square foot, 4% square foot, and 46 square foot in the same 
field 

1-SQUARE-FOOT UNITS 


Average wireworms obtained from samples of the specified number of 








units 
or a 
Field no. a * = 
Relative Relative Relative 
Mean standard Mean ! standard Mean standard 
error error error 

Number | Percent | Number Percent Number Percent 
1 0. 22+0). 04 18. 1 
2 0. 283+0. 09 32.1 | 0.3240. 08 25.0 09 20.0 
68+ .17 25.0 78+ .18 | 23. 1 11 15.9 
j 1. 08+ x0 21.8 _ B+ 18 20.5 11 | 14.9 
92+ .20 21.7 96+ .16 16.7 12 14.8 
6 444 14 31.8 824 .17 20.7 ‘ 4 14. 6 
7 1. 22+ .35 23.7 | 1.244 .22 17.7 | 1.4 18 14.4 
8 1.484 .31 20.9 | 1.444 .25 17.4 1.7 2) 11.4 
9 2.56 . 19.5 | 3.144 42 13.4 | 2 24 8.1 
10. 3.322 . 51 15.4 | 2.964 .36 12.2 | 2.4 34 11.4 
11 3. 24+ .52 16.0 | 3.10+ .38 12.3 | 3 33 10.4 
12 3.004 .72 24.0 | 3.184 44 13.8 | 3. 30 9.4 
13 100+ .79 19.8 | 3.66+ .49 13.4 | 3.924 .39 9.9 
4 4. 24+1. 00 23.6 | 3.884 .43 t1. 8 3.944 .38 9.6 
15 3.202 .4 13.8 | 3.404 .37 10.9 | 3.694 .32 8.7 
16 4.124 .48 11.7 | 4.024 .40 9.9 | 4.254 .20 4.7 
17 4.00+ . 51 12.8 | 4.364 .43 9.9 | 471+ .37 7.9 
Ig 6. 96+1. 09 15.7 | 5. 66+ . 66 11.7 6. 014 47 1.8 
19 5. 80+ .70 12.1 | 6.76 .80 11.8 | 6.314 . 52 &.2 
29 6.96+ .90 12.9 | 7.084 .93 13.1 | 7.02> . 59 8.4 
21 8. 5641. 30 15.2 | 7.944 . 8 10.1 | 7.534 . 56 7.4 
22 9. 28+1. 19 12.8 | 8.964 .71 7.9 | 7.964 .51 6.4 
23 8. 6041. 23 14.3 | 7.864 . 56 7.1 | 8.204 .47 5.7 
24 11. 121. 63 14.7 | 7.604 .77 10.1 | 8. 83+ .71 8.0 
25 14. 1622. 12 15.0 | 9.46+ .84 8.9 |10.78+ .76 | 7.1 

44-SQUARE-FOOT UNITS 
0. 20-40. 10 50.0 | 0. 1240.05 41.7 | 0.12+0.04 33.3 
2 16+ .07 44.8 -124 .05 41.7 12+ .03 25.0 
3 40+ .10 25.0 . 20+ .06 30.0 29+ .05 17.2 
4 60+ .15 25. 0 52+ .10 19, 2 47+ .0% 14.9 
5 Sit .12 21.4 
6 1.40+ .29 20.7 | 1.004 .16 16.0 77+ .11 14.3 
r) 108+ . 26 24.1 | 1.144 .19 |} 16.7 90+ .11 12.2 
& 1.68+ . 37 22.0 | 2.08+ . 26 12.5 | 1.864 .19 10. 2 
Kise SQUARE-FOOT UNITS 
| 0. 120. 06 50.0 | 0. O80. 04 50.0 | 0.0540. 02 40.0 
2 OR+ .04 50.0 10+ .03 30.0 
5 12+ .08 66.7 - 10+ .05 50.0 -ll+ .08 27.3 
4 16+ .07 43.7 18+ .06 33.3 25+ .06 24.0 
5 644% .20 31.3 -464 .10 21.7 41+ .06 14.6 
6 52+ .13 25.0 74+ .18 24.3 50+ .10 | 20.0 
7 60+ . 16 26. 6 624 .13 21.0 64+ .09 14. 1 
8 72+ .16 22. 2 . 664 .11 16.7 73 .09 12.3 
9 764 .15 19.7 84+ .12 14.3 88+ .10 11.4 
10__. . — . 1.00+ . 24 24.0 
_s™= 1.04% . 29 27.9 | 1.084 .22 20.4 | 1.134% .14 12.4 
12.. “ pis ee Sees es . “ 2.004 .31 15.5 
13_. Shnesuntieremnenheee 1.802 .41 | 22. 8 2.144 . 28 | 13.1 2.19% .24 11.0 
| | j 





See footnote table 2. 
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To determine the smallest number of units that may be taken to 
give a mean with a relative standard error considered reasonable for 
wireworm populations, comparisons were made on the basis of four 
different population densities for each of the three sampling units. 
Means and standard deviations were calculated from 100 units of 
each size from the four different population densities. Considering 
the standard deviations obtained from these distributions as repre- 
sentative of the true standard deviations of the population, theoretical 
standard errors of the mean were computed by varying the number 
of units in the sample from 1 to 100. These errors were converted 
into relative errors of the mean to facilitate comparison and then 
plotted in figure 7. 

A study of these graphs reveals the more important characteristics 
shown by the data. Little practical effect upon the sampling error 
is obtained at any population level by using mere than 50 units to 
determine the average number of wireworms per square foot of soil, 
but at least 50 units are necessary at low population levels. There 
is demonstrated also a decreasing necessity for taking 50 units, as 
the size of the sampling unit is increased from \. square foot to 1 square 
foot. At the intermediate and high levels of population a smaller 
number of units is sufficient. 

Except for low population levels, the number of units required for 
a 20- to 40-percent error range is not over 25 for 1-square-foot and 
\4-square-foot units, or 50 for \.-square-foot units. 


DISCUSSION 


There are certain limits to the amount of sampling that is practicable 
in entomological field surveys. The field man has only a limited 
amount of time to devote to sampling. The work here reported 
shows that fifty 1-square-foot units, sifted, are enough to give fairly 
accurate estimates of the number of wireworms at any population 
level. These 1-square-foot siftings, taken annually in the same field, 
permit comparisons of wireworm populations from one season to 
another. The degree of crop damage done by different numbers of 
wireworms can be determined. The 1-square-foot unit is large 
enough to be used where depth records of the wireworms are desired 
in special studies of vertical migration. For purposes pertaining to 
the life history, however, where counts of eggs and young larvae are 
required, or for determining the size of wireworm infestations during 
summer control practices, the }.s-square-foot units, washed, are 
necessary. Fifty of these units, washed, will give fairly accurate 
estimates of the population, except for very low population levels, 
where 100 such units may be necessary. 

Each of these sampling units has certain advantages. The tech- 
nique developed for obtaining wireworm counts by sifting calls for 
more costly equipment than that for obtaining counts by washing. 
With each method the labor of two men is required. From 30 to 50 
cubic feet of soil per day can be sifted with the power sifter, whereas 
one hundred \¢-square-foot units can be washed and examined in a 
working day. 

SUMMARY 


Practical quantitative methods of estimating wireworm populations 
in soil have been presented. Efforts have been made to find a time- 
and labor-saving system of sampling that will give fairly accurate 
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results. In brief, the methods selected consist of taking small sam- 
pling units of soil at random over the field. The field equipment 
consists of a metal sampling form, a soil-sifting machine or a soil- 
washing outfit, a tamper, a square-edged shovel, 10-quart pails, and 
record blanks. 

The data obtained with sampling units of three sizes—1l, 4, and 
Ne square foot—have been subjected to a statistical analysis. The 
wireworm distribution is moderately asymmetrical for the 1-square- 
foot unit, the class frequency decreasing with markedly greater 
rapidity on one side of the maximum than on the other. The distri- 
butions show close agreement with the Poisson series. It is shown 
that the 1l-square-foot units give more precision in estimating the 
mean than either of the smaller units. 

In determining the smallest number of sampling units that may be 
taken, comparisons were made on the basis of four different population 
densities for each of the three units. Little practical reduction of 
the sampling error is obtained at any population by using more than 
50 samples to determine the average number of wireworms per 
square foot of soil. At low population levels it is necessary to take 
at least 50 samples. At the intermediate and high levels of population 
a smaller number of samples is sufficient. 
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THE BIOLOGICAL VALUE OF THE PROTEINS IN HEGARI 
AND THE SUPPLEMENTAL VALUE OF CERTAIN PRO- 
TEIN CONCENTRATES USED IN FARM ANIMAL FEED- 
ING' 


By MARGARET CAMMACK Situ, nutrition chemist, and GLApys HartTLey RoeumM, 
assistant nutrition chemist, Arizona Agricultural Experiment Station 


INTRODUCTION 


Previous work in this laboratory has shown that the sorghum grain 
hegari, like other grains, is not adequate as the sole source of protein 
in the diet of rats (14).2 That the protein in grains must be supple- 
mented has long been recognized, and it is common practice to incor- 
porate in farm animal rations made up largely of grains, a protein- 
rich food such as tankage, cottonseed meal, etc. 

Osborne and Mendel (13) and McCollum, Simmonds, and Parsons 
(1, 2, 3, 4, 5) have demonstrated the way in which proteins supple- 
ment each other in nutrition, and have pointed out that the nutritive 
value of a protein when fed as the only source of protein in the ration 
is not a measure of its nutritive value when fed in combination with 
other proteins. It is a matter of great practical importance, therefore, 
so to combine natural foods as to secure protein mixtures of high bio- 
logical value. Thus the question of the best protein concentrate from 
the standpoint of supplemental value to hegari proteins arises. This 
paper gives the results of an investigation of the nutritive value of 
hegari proteins and the extent to which certain protein concentrates 
exert a supplemental effect. 


IXPERIMENTAL PROCEDURE AND RESULTS 
GROWTH EXPERIMENTS 


In the first series of experiments the nutritive value of hegari as 
the sole source of protein in a ration, and the nutritive value of hegari 
supplemented by certain protein concentrates or amino acids, was 
studied by means of growth experiments. The method followed in 
part of this investigation was essentially that used by McCollum, 
Simmonds, and Parsons (1, 2, 3, 4,5). Albino rats were taken at the 
weaning age (4 weeks) when weighing 45 to 65 g, placed in separate 
cages, and given a ration containing a single protein food or a mix- 
ture of proteins at a level of 9 percent, for McCollum considered this 
a critical level even when the protein was of excellent quality. The 
diets were so constituted (table 1) that all the other known dietary 
essentials were supplied in ample amounts, so that the growth re- 
sponses would be a measure of the nutritive value of the protein frac- 
tion which was the only variable component. Apparently all the 
rations were palatable to the rats except those containing tankage or 
meat scrap. The food intake on these rations was low, which may 
have influenced the utilization of the proteins for growth to some 
extent. 

"1 Received for publication June 30, 1936; issued March 1937. 


? Reference is made by number (italic) to Literature Cited, p. 145. 
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The basal diet and distilled water were available to the animals at 
all times, and records were made weekly of the food consumed and of 
the body weights for an 8-week test period. 

The protein-rich materials whose value as supplements for hegari 
proteins was studied included cottonseed meal, soybean-oil meal,* 
dried skim milk, tankage, gelatin, wheat gluten, dried egg albumen, 
and alfalfa-leaf meal. In each case the protein mixture was fed at a 
9-percent level, two-thirds of whch came from hegari and one-third 
from one of the above-mentioned high-protein foods (table 1). 

Table 2 summarizes the results obtained in terms of the average 
gain in weight and average total food consumed in the 8-week period, 
and the average gain per gram of protein eaten by the rats on each of 
the experimental rations. Great differences in value for the support 
of growth of the various protein foods when fed singly at this critical 
protein level may be seen. The best gains in weight per gram of food 
or protein consumed were obtained from dried skim milk and from 
egg albumen. Cottonseed meal and soybean-oil-meal proteins were 
somewhat superior to the proteins of wheat gluten, alfalfa-leaf meal, 
gelatin, and tankage for the support of growth. 

Distinct supplementary relations between certain of these proteins 
and hegari proteins were also observed when they replaced part of the 
hegari proteins in the ration (i. e., 6 percent of hegari protein and 
3 percent of dried skim-milk protein 1). A comparison of experi- 
mental gains in weight per gram of protein consumed and the calcu- 
lated gains in weight which would be expected from a mixture of 
hegari and the protein concentrate in the proportion fed, revealed a 
definite supplemental relation in some cases. For example, mixtures 
of hegari with dried skim milk, or gelatin, or egg albumen, or tankage, 
or alfalfa-leaf meal provide protein more efficient for the promotion 
of growth than could be expected if no supplementary relationship 
existed. It will be seen that when hegari is mixed with dried skim 
milk, gelatin, egg albumen, alfalfa-leaf meal, or tankage, the actually 
obtained gains per gram of protein mixture are greater than the 
weighted mean gain per gram of protein mixture calculated from the 
gain per gram of each protein in the mixture when fed singly. It 
would appear therefore that there is a portion of the protein fraction 
of hegari or of the protein concentrate which is unused when fed 
singly, but which becomes available when fed in conjunction with 
each other. The nutritive value of the proteins of gelatin is very poor 
when fed as the only source of protein in the ration, so much so that 
rats on this ration cannot even maintain their weight. However when 
gelatin is combined with hegari, the growth-promoting value of the 
protein mixture is far superior to either protein when fed singly. The 
supplementary relations are practically as great between gelatin and 
hegari proteins as between dried skim milk and hegari proteins. 
However, it is true, of course, that because of the extremely poor 
quality of gelatin proteins, the resulting mixture of gelatin and hegari 
protein at the 9-percent level is markedly inferior in nutritive value to 
the hegari-milk protein mixture. 


The soybean-oil meal used in these tests was a well-cooked, hydraulic-process product. 
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TABLE 2.—Comparative efficiency for the support of growth of hegari proteins, when 
fed at 9 percent protein level for 8 weeks, alone and when supplemented by certain 
protein concentrates 


Average total Calcu- 
intake lated 
weighted 


——_, | Observed 


| Aver- | Aver- | | gain in Pa ny 
: : . age age weight tee he 
5 2 pin (¢ nate values) ats tee wag t 
ource of protein (approximat alues Rats initial | gain in per gram b.. 3 
weight | weight Food Pro- | of protein or gram 
} tein ingested | POT &te 





of protein 
mixture 
ingested ! 
Number, Grams | Grams | Grams | Grams | Grams Grams 
9 percent from hegari 10 55 24 452 40. 68 +0). 59 
9 percent from cottonseed meal —_ 4 61 80 | 550 49. 50 +1. 62 
6 percent from hegari, 3 percent from 
cottonseed meal... . wane 1 56 69 657 59. 13 +1.17 +0). 93 
9 percent from dried skim milk_. 4 52 162 677 60. 93 +2. 65 
6 percent from hegari, 3 percent from 
dried skim milk ms 4 52 174 792 71. 28 +2. 44 +1. 28 
9 percent from aifalfa- leaf meal 2 58 18 407 36. 63 +. 49 
6 percent from hegari, 3 percent from 
alfalfa-leaf meal nese 4 60 85 655 58. 93 +1. 44 +. 56 
9 percent from soybean-oil meal : 4 59 40 377 33.9 +1.17 
6 percent from hegari, 3 percent from 
soybean-oil meal 4 60 70 540 48. 6 +1. 44 +.78 
9 percent from wheat gluten = 2 58 50 488 43. 92 +1.14 
6 percent from hegari, 3 percent from 
wheat gluten F 2 57 28 350 31. 50 +. 88 +.77 
9 percent from egg albumen 4 63 8S 406 | 36.54 +2. 40 
6 percent from hegari, 3 percent from e a4 
albumen 4 56 135 700 63. 00 +2.14 +1.19 
9 percent from gelatin 4 52 20 248 22. 32 +. 89 
6 percent from hegari, 3 percent from 
gelatin 4 56 46 466 41.94 +1. 10 +. 09 ‘ 
6 percent from hegari, 3 percent from i 
tankage 4 55 7 190 | 17.1 +. 41 
® percent from tanks age 4 58 16 139 2. 53 +1. 27 +. 53 ; 


1 Expected on the assumption that there is no supplemental effect among the proteins in the mixture 


On the other hand, no appreciable supplementary relation between 
the proteins of hegari and wheat gluten or cottonseed meal could be 
observed. The value of the incorporation of these proteins in hegari 
rations would depend only upon an increase in the protein content or 
quality but not upon any supplementary relationship. 

By using the perhaps more quantitative paired-feeding procedure 
considered desirable by Mitchell, information as to the supplemental 
effect of certain amino acids has been obtained. For this study the 
rats were carefully matched as to sex, litter, and initial weight and the 
food intake of one member of a pair was restricted to the food intake 
of the other. The rats were weighed weekly and accurate records of 
the food intake were made daily for a period of 6 weeks. Using the 
comparative growth response in from 7 to 10 pairs of rats as the cri- 
terion, the supplementary relations between hegari and cystine, trypto- 
phane, and lysine were investigated. The amino acid additions in the 
amount of 0.02 g fed daily were made to the hegari ration used in the 
previously described experiments. (See table 1 

The results are presented in table 3, for tryptophane, cystine, and 
lysine, respectively, and the effect of feeding gelatin as a protein high 
in lysine is also shown. In the last case, the gelatin provided 3 per- 
cent and the hegari 6 percent of the protein in ‘the ration. Inspection 
of these tabulated results shows that no supplemental effect between 
cystine and hegari proteins or between tryptophane and hegari pro- 
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teins could be consistently demonstrated. However, in every case, 
the addition of lysine accelerated the growth rate of the rats receiving 
this supplement to their hegari ration. Again, it may be noted in 
table 3 that when gelatin as a source of lysine was incorporated in the 
i ration, the growth of the gelatin-fed member of the rat pair exceeded 
the growth rate of its pair mate in every case, thus giving further evi- 
dence that lysine is a limiting factor in the hegari-protein make-up. 


TABLE 3.—-Effect of the addition of tryptophane, cystine, lysine, and gelatin to hegari 





proteins, as shown by paired-feeding growth experiments of 6 weeks’ duration 
TRYPTOPHANE 
Initial | Final | Total oy Supple- 
Pair no. Source of protein weight | weight food of ms mental 
of rats | of rats | intake m™.o 1, effect 
Grams | Grams | Grams | Grams 
{ Hegari ee ace ‘ 47 264 3 P 
\Hegari +0.02 g tryptophane daily - -- 42 43 263 l ~ 
. { Hegari aig rex, aS 60 69 295 9 ns 
. ~|\Hegari +0.02 g tryptophane daily. -- e 48 56 298 8 _ 
\f Hegari- ... eis 48 49 254 1 
\Hegari +0.02 g tryptophane daily 48 52 253 4 . 
P { Hegari-_. ee ie 3 53 58 224 5 |. ce 
\Hegari +0.02 g tryptophane daily--_-__- 55 61 224 6 + 
{ Hegari om - > a Re ahs 60 69 267 | Ce 
\Hegari +0.02 g tryptophane daily. ---- 60 69 265 9 0 
, { Hegari ent 3 ‘: ces 63 74 271 _. =e 
: \Hegari +0.02 g tryptophane daily. .__.-. 61 75 271 14 + 
» { Hegari " = 55 66 258 ll 
. i — \Hegari +0.02 g tryptophane daily - _- 68 76 256 8 _ 
8 { Hegari te 9 71 270 12 ‘ i 
\Hegari +0.02 g tryptophane daily 58 67 270 ) 
{ legari a 7.7 ; 
Average \Hegari +0.02 g tryptophane daily cd 7.4 in 
CYSTINE 
. f Hegari-_--_-.--- ae ey 63 73 268 10 aa 
; ' \Hegari +0.02 g cystine daily 61 70 268 9 _ 
. { Hegari : x 57 68 265 |. : 
. \Hegari +0.02 g cystine daily 56 71 264 15 + 
> j Hegari ‘ 70 81 254 11 wits 
‘ \Hegari +0.02 g cystine daily 69 77 258 8 vt 
F { Hegari 65 74 267 yg ~ 
\Hegari +0.02 g cystine daily 65 75 264 10 4 
{ Hegari 59 67 260 8 , 
\Hegari +0.02 g cystine daily 58 68 253 10 + 
, { Hegari 59 76 282 17 
; \Hegari +0.02 g cystine daily ; 60 75 284 15 _ 
° { Hegari : 60 77 277 17 |.. 
: \Hegari +0.02 g cystine daily | 59 71 277 12 - 
{ Hegari ian . ‘ es hee 11.9 
Average | Hegari +0.02 g cystine daily | ee ae. 
LYSINE 
{ Hegari 60 45 187 —15 
\Hegari +0.02 g lysine daily- 60 58 186 —2 + 
. j Hegari 56 62 229 6 
— \Hegari +0.02 g lysine daily 57 71 232 14 - 
. { Hegari 55 59 257 4 
\Hegari +0.02 g lysine daily- 55 65 262 10 . 
4 { Hegari 60 62 241 2 
— \ Hegari +0.02 g lysine daily 61 73 240 12 + 
i { Hegari * —.8 
Average \Hegari +0.02 g lysine daily . 8.5 + 





! Minus sign (—) indicates a loss. 
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TABLE 3.—Effect of the addition of tryptophane, cystine, lysine, and gelatin to hegari 
proteins, as shown by paired-feeding growth experiments of 6 weeks’ duration—Con. 


GELATIN 





_ - 4 G 
; | Initial | Final | Total or ol | Supple- 
Pair no. Source of protein weight | weight food Ps mental 
of rats | of rats | intake weight effect 


Grams | Grams | Grams Grams 
283. 5 14 


7 f Hegari 55 69 
\Hegari +0.02 g gelatin daily- 56 78 281.0 22 + 
° j Hegari _ - 55 74 310.0 19 
™ \Hegari +0.02 g gelatin daily 53 78 309.0 25 + 
3 { Hegari 60 80 313.0 20 
. oe +0.02 g gelatin daily__- 61 84 315.0 2 + 
{ jHeg: 48 | 64 | 268.0 16 | 
\Heeari +0.02 g gelatin daily 5O 75 266. 5 25 + 
{ Hegari 55 65 772.0 10 
\Hegari +0.02 g gelatin daily 53 70 274.0 17 + 
6 { Hegari 58 68 284.0 10 
: \ Hegari +0.02 g gelatin daily 58 71 281.0 13 | + 
- { Hegari 60- 65 254.0 | 5 : 
¥ \Hegari +0.02 g gelatin daily 59 | 72 | 252.5 13 + 
8 { Hegari 61 73 276.0 12 
\Hegari +0.02 g gelatin daily- 60 80 278.0 20 + 
9 { Hegari : 57 71 307.0 14 
; \Hegari +0.02 g gelatin daily 56 76 | 301.0 20 | + 
10 | Hegari 59 75 | 311.0 16 
\Hegari +0.02 g gelatin daily 59 78 307.0 19 + 
Average { Hegari 13. 6 


\Hegari +0.02 g gelatin daily 19.7 + 


NITROGEN-BALANCE EXPERIMENTS 


In the second series of experiments planned to determine the nutri- 
tive value of the proteins of hegari and its supplementary relations 
with certain protein-rich foods, the more direct method of investiga- 
tion by means of nitrogen-metabolism studies was followed. 

This method, based upon that used by Karl Thomas, has been 
worked out in great detail by Mitchell (6, 8) for use with rats. It 
consists essentially in measuring the percentage of the absorbed food 
nitrogen that is retained in the body, the result being referred to as 
the biological value of the protein. 

Male albino rats only were used in the experiments. They were 
housed in galvanized-iron cages with raised screen bottoms placed 
upon porcelain plates. Distilled water was provided in a glass bottle 
attached to the outside of the cage. The analyzed food rations were 
placed in nonscatter glass cups and a quantitative record of the food 
intake was made. The urine was absorbed on large circles of filter 
paper which had previously been soaked in a 1 to 1 solution of hydro- 
chloric acid and dried. The filter paper with absorbed urine was 
changed daily and the feces collected quantitatively and preserved 
with concentrated sulphuric acid. At the end of the 5-day test 
period, which had been preceded by a 3-day fore period for dietary 
adjustment, the filter papers were washed free from urine with hot, 
dilute hydrochloric acid, made up to volume, and an aliquot was 
taken for nitrogen analysis by the Kjeldahl method. All the feces 
obtained in the test period were digested with concentrated sulphuric 
acid, mercury, and potassium sulphate, diluted and made up to volume, 
and an aliquot was analyzed for nitrogen 

For each rat, the endogenous nitrogen in urine and in feces (per 100 g 
of body weight and per gram of food consumed respectively) was 
determined on the low-nitrogen whole-dried egg diet recommended by 
Mitchell and Carman (9), except that the lard was omitted and the 
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percentage of carbohydrate correspondingly increased. As _ hegari 
rations supplying approximately 9 percent of hegari proteins must of 
necessity be largely hegari, this adjustment was made in order that 
the energy value of all the experimental rations used for nitrogen- 


balance studies should be approximately the same. 
The experimental rations were again those shown in table 1. They 
: were so planned as to be approximately the same not only in energy 


value but also in roughage content. 

The results of these nitrogen-balance studies expressed as biological 
value of the proteins when fed singly and also in mixture with hegari 
proteins appear in tables 4 and 5. 

It may be seen that hegari proteins have the relatively low biological 
value of 53.2 percent. About half of the absorbed protein in hegari, 
therefore, is not utilizable for growth, because the assortment of amino 
acids essential for synthesis into the complex body tissue is not com- 
plete. Thus hegari protein is capable of being supplemented to a large 
extent by those proteins that supply the amino acids which are de- 
ficient in the hegari protein make-up. 

The biological value of dried skim milk, alfalfa-leaf meal, cottonseed 
meal, and soybean-oil meal as determined in this laboratory agree 
closely with the results of similar measurements made by Mitchell and 
others (7, 10, 11, 13, 15). Determination of the biological value of 
these protein foods when fed singly and in combination to the same 
group of animals leaves no doubt as to the reliability of the observed 
supplemental values. 

The supplementary relations between hegari protein and the pro- 
teins of dried skim milk, cottonseed meal, meat scrap, alfalfa-leaf 
meal, and soybean-oil meal are shown in table 5. The differences 
between the true or observed biological value of the mixture and the 
calculated biological value, which is the weighted mean of the biological 
values of the proteins in the mixture when fed separately, gives the 
supplemental value of the proteins concerned. 





TaBLE 4.—Biological value of proteins from various sources when fed singly and in 
mixture with hegari proteins 


HEGARI PROTEINS FED AT A 9 PERCENT LEVEL 
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12406 162 165 | 14.0 | 180.8 | 26.94 7. 90 173.0 | 105.48 | 79.9 93. 1 53. 8 
12408 163 164 | 15.2 | 196.3 | 38.46 | 11.26 185.0 | 112.95 | 0.3 104.7 56. 6 
12409 156 159 | 18.0 | 232.3 | 44.99 | 32.07 200.0 | 132.66 | 95.8 104. 2 52. 1 
12410 149 149 | 13.6 | 175.6 | 27.13 6.0 169.6 | 100.46 | 68.7 100.9 59. 5 
11893 197 201 14.8 | 190.9 | 48.4 17.8 173. 1 112.6 94. 1 79.0 | 45. 6 
12246. 235 238 | 13.4 | 172.9 | 34.4 11,9 161.0 | 102.5 62.8 98. 2 60.9 
12415 254 261 17.0 | 219.3 | 45.4 19.4 199.9 | 128.2 85.3 114.6 57.3 
12418 212 217 15.6 | 201.2 | 40.9 15.0 186.2 | 129.9 92.4 93.8 50.4 
8659 151 159 | 10.2 | 142.6 | 36.5 30.6 112.0 91.3 55.1 56.9 50.8 
8676 158 159 7.8 | 109.0 | 25.2 17.6 91.5 75.5 40.0 51.5 56.3 
8678 159 160 8.0 | 111.8 | 26.8 19.4 92.4 74.3 38. 0 54.5 59. 0 
S684 147 152 8.8 123.0 | 28.1 21.8 101.2 79.5 46.9 54.3 53.7 
S6R5_ _. 150 151 6.8 95.1 | 23.9 20.5 74.6 63. 1 30. 1 44.5 59.7 
12111 . 282 285 | 11.0 | 153.8 | 37.6 27.4 126.4 | 1901.5 62.7 63.7 51.4 
12112 276 275 | 19.2 | 128.6 | 36.2 23.8 104.8 93. 2 55.7 49.1 46.9 
12113 261 259 8.8 | 123.0 | 31.8 20.5 102. 5 93.0 55.5 47.0 45.9 
12115 255 252 9.8 | 139.1 | 37.6 26.3 112.8 96. 6 58. 6 54. 2 48.0 
12114... 292 | 292 | 10.8 | 151.0 | 45.0 | 32.7 112.8 | 98.8 | 60.6 57.8 | 48.9 
Average ; 53.2 





142 


Journal of Agricultural Research 





Vol 
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[aBLE 4.—Biological value of proteins from various sources when fed singly and in 
mixture with hegart proteins—Continued 


ALFALFA-LEAF-MEAL PROTEINS FED AT A 9PERCENT LEVEL 
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12408 ; é 209.0 |209.0 | 11.0 | 158.4 20. 2 ant Mee 85.7 | 50.2 108.2) 68.3 
2410 180.0 |181.0 9.2 | 132.5 56. 4 | 43. 1 89.4 91.3 | 34.6 54.8 61.3 
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12411 210.0 |208.0 7 51.5 31.6 167. 1 95.8 | 53.8 113.3 67.8 
11983 200.0 |200.0 2 40.8 14.3 164.9 99.1 | 66.1 | O88 59.9 
8676 181.0 |192.0 4 22.7 15, 2 169. 3 99.9 | 64.4 104.9 62.0 
A verage_..- — — otuwe bebe - ae 61.2 
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2406 225. 0 19.6 | 245.0 53. 4 36.8 208. 2 111.4 77.9 130.3 62.6 
: 12407. . - 225. 0 22.2 | 277.5 59. 2 6.9 270.6 132.2 | 91.3 179.3 66.3 
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; 12409. “ -/212.0 22.2 | 277.5 57.9 42.6 234.9 130.5 | 80.7 154. 2 65.6 
L 12410 180. 0 16.6 | 207.5 57. ¢ 33. 2 174.3 100.2 | 43.5 130.8 75.0 
2411... 1174.0 17.4 | 217.5 41.3 19.9 197.6 101.4 | 58.6 139. 0 70.3 
12408 al 248.0 18.4 | 231.8 70.6) 35.5 196. 3 101.9 | 59.7 136. 6 69. 6 
t 12411 . 205. 0 15.0 | 189.0 | 56.0 37.5 151.5 91.8 | 41.4 110.1 72.7 
12421... eesti 200.0 11.2 | 156.8 29.7 13.0 143.8 | 82.9 | 49.7 94.1 65.4 
S684... 170.0 |166.0 9.5 | 147.2 21.6 15.2 132.0 | 80. 1 47.6 84.5 64.0 
Average... , . sess . : pawn lonneon -| 67.7 
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12411 2.0 |191.0 9.8 | 185.2 39.8 159. 5 115.9 | 86.1 73.4 46.0 
12407... 9.0 (241.0 | 15.8 | 208.6 90.0 245.9 174. 2 |128.9 117.0 47.5 
12408 _.._. 0 (230.0 | 16.6 | 313.7 91.4 254. 0 182.0 |142.4 111 6 43.9 
12409_. 5.0 |235.0 | 16.8 | 317.5 90. 2 238.9 182.6 |127.4 |) 1115 46.7 
14206 3.0 |233.0 | 13.8 | 260.8 59.1 222.4 | 152 117.1 105.3 | 47.3 
11983 3.0 |198.0 12.0 | 226.8 66.8 185.1 | 143.5 |109.0 76.1) 41.1 
12407 9.0 |264.0 12.2 | 230.6 58.8 194.0 151.8 |103.7 90.3 46.5 
12246 3.0 |231.0 | 12.0 | 226.8 41.1 206. 3 152.9 [113.5 92.8 | 45.0 
12415 255.0 |235.0 7.6 | 143.6 26. 2 129.0 118.2 | 75.1 53.9 | 41.8 
12418 9.0 203.0 | 11.8 | 223.3 43.6 198. 7 149.0 |112.0 86.7 43.6 
12421_. .0 |201.0 | 11.4 | 215.5 37.3 195, 2 133.8 |100.3 | 95. 0 | 48.7 
8659... 38.0 |157.0 | 14.0 | 264.6 | 70.8 219.8 | 159.0 112.7 | 107.1 | 48.7 
8676... 31.0 160.0} 82 | 121.6| 33.8 95.4| 804/448] 506] 53.0 
3678... 39.0 |161.0| 7.2| 106.8] 36.6 | 77.7 | 72.6|362| 41.4| 533 
S684 56.0 |152.0 | 8.0| 118.6] 34.8 | 9.1} 86.1 43.5) 46.6| 51.7 
$085 ae 0 |148.0| 5.0] 74.2] 227 56.2| 59.6/ 26.6] 20.6| 52.7 
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TABLE 4. 


mixture with hegari proteins—-Continued 
MIXTURE OF HEGARI AND MEAT-SCRAP PROTEINS FED AT 


34 HEGARI PROTEIN AND 146 MEAT-SCRAP PROTEIN 
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12406 238.0 |247.0 | 18.6 | 186.0 3. : 27.5) 158.5 | 111.4! 75.9 82.6} 52.1 
12407 250.0 256.0 | 20.2 | 202.0) 36 0 | 202.0 107.5 | 62.0 140.0 | 69.3 
612408 230.0 |238.0 | 20.0 | 200.0 5. 5 8.3 | 191.7 | 103.0 | 63.9 27.8 | 66.7 
11983 2: 200.0 200.0 | 15.8 | 158.0 47.0 14.3 | 143.7 96.9 | 62.9 80.8 | 56.2 
12246 235.0 |240.0 | 16.6 | 166.0) 50.7 | 22.1] 143.9) 100.8 | 61.1 82.8 | 57.5 
12418 213.0 |216.0 | 17.4 | 174.0 56. 6 25. 8 148.2 | 108.6 | 70.9 77.3 | 52.2 
12421 204.0 209.0 | 16.2 | 162.0 ) 52 26.0 136.0 | 82.3 | 48.4 87.6 | 64.4 
SO6N4 158.0 (161.0 | 10.0 | 139.5 28. : 22.0 117.6 75.3 | 42.7 74.9 | 63.7 
8659 165.0 |164.0 | 8.2 | 114.4) 32. 26. 2 88. 2 70.8 | 34.6 53.6 | 60.8 
8676 167.0 |167.0 | 7.6 | 106.0] 329] 24.3 81.7 64.2 | 28.6 53. 1 65.0 
Average a 60.8 

SOY BEAN-OIL-MEAL PROTEINS FED AT A 9-PERCENT LEVEL 

15008 152.0 |157.0 7.2 | 109.5 21.7 103.7 80.9 | 37.2 66.5 | 64.1 
16088 116.0 |122.0 | 10.9 | 166.8 | 33.2) 21 145.4 60.4 | 49.6 95.7 | 65.8 
16089 118.0 [127.0 | 11.0 | 168.3 | 32.6) 21 146.8 60.6 | 47.9 98.8 | 65.2 
15007 192.0 |200.0 | 12.7 | 194.3 | 15.2 | 193.8 84.3 | 57.8 136.0 | 70.2 
16191 113.0 |117.0 | 10.2 | 154.5 | 26.5 18 138.7 64.6 | 48.2 90.5 | 65.2 
16192 86.0 | 94.0) 6.0) 91.8) 9.58 91.8 64.2 | 41.2 50.6 | 55.2 
Average _. nis 64.3 
MIXTURE OF HEGARI AND SOYBEAN-OIL-MEAL PROTEINS FED AT A 9PERCENT 


LEVEL, 74 HEGARI PROTEIN AND 








15008 166.0 |170.0 | 18.0 | 278.0 


50.6 | 10.7 | 267.3 | 128.8 
16088 111.0 |111.0 | 10.3 | 175.0 32.7 21.6 | 153.4 163. 9 
16089 105.0 109.0 | 11.0) 170.0} 30.7 19.6 | 150.4 76.8 
16044 82.0 | 86.0| 6.3 | 97.0 17.3 10.2 | 86.85 | 55.95 
15007 192.0 |192.0 | 12.9 | 193.3 | 33.1 12.2 | 181.1 92.9 
16091 _. 125.0 |134.0 | 13.4 | 207.0] 29.6 2.9 | 204.1 66.6 
16191 98.0 |103.0 | 12.4 | 191.6 28.3 17.8 | 173.8 | 59.7 
16192 | 82.0 | 94.0 | 10.6 | 163.8 22.6 0.7 | 163.1 | 59.7 


Average 


TABLE 5. 
meal, soybean-oil meal, dried skim milk, meat scrap, 
by nitrogen-balance experiments 
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from 
from 
from 
from 
from 
from 
from 


cottonseed meal _ ad 

hegari, 3 percent from cottonseed meal 
dried skim milk__- 

hegari, 3 percent from dried skim milk 
alfalfa-leaf meal_-_... tos 7 
hegari, 3 percent from alfalfa-leaf meal 
soybean-oil meal ss iaaicbatoes 
6 percent from hegari, 3 percent from soybean-oil meal can 
9 percent from meat scrap . we 
6 percent from hegari, 3 percent from meat scrap. - - 


¥ percent 
6 percent 
9 percent 
6 percent 
9 percent 
6 percent 
9 percent 








1+ Probable error 


4 SOYBEAN-OIL-MEAL PROTEIN 


81.2 | 186.1] 69.6 
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67.2} 113.9] 62.9 
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} 41.6) 121.4] 74.5 
= 67.9 


Supplementary relations between proteins of hegari and those of cottonseed 
or alfalfa-leaf meal as shown 
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56. 6 0.6 
63.3 14.0 
"55.9 11. § 
57.0 10.9 
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Thus it is evident that dried skim milk, alfalfa-leaf meal, soybean- 
oil meal, and meat scrap exhibit supplementary value in the order 
given, whereas cottonseed meal has no appreciable supplemental 
effect. The use of cottonseed meal as a protein concentrate in rations 
; made up largely of hegari serves therefore only to increase the total 
. protein in the diet. 

Dried skim milk, alfalfa-leaf meal, and meat scrap, on the other 
hand, may be used not only to increase the total protein content of 
the ration but also to make available an otherwise nonutilizable por- 
tion of the proteins in the mixture, thus serving in two ways to in- 
crease the protein in the ration. It would appear that the most effi- 
cient proteins from the standpoint of supplemental value for hegari 
proteins are those rich in lysine. 






CONCLUSIONS 


The following conclusions have been drawn from growth- and nitro- 
cen-balance experiments with rats fed rations in which hegari served 
as the only source of protein in an otherwise adequate ration, or was 
partly substituted with certain protein concentrates: 

(1) Hegari proteins are inadequate for the support of growth, as 
shown by a relatively low rate of gain per gram of protein consumed 
by the rats used in the growth studies, and a low biological value of 
the protein as determined by nitrogen-balance experiments. 

(2) The protein-rich foods which have shown definite supplemental 
relations with hegari proteins include dried skim milk, alfalfa-leaf 
meal, soybean-oil meal, meat scrap, gelatin, and tankage. 

(3) The protein concentrates studied which were found to be with- 
out supplemental relations were cottonseed meal and wheat gluten. 

(4) Paired feeding experiments showed that the amino acid, lysine, 
has a marked ability to supplement hegari proteins, whereas trypto- 
phane and cystine have not. It would appear that lysine deficiency 
is at least partly responsible for the low biological value of the pro- 
teins of hegari and that the proteins which are rich in this amino acid 
have a strong supplemental value for hegar. 
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PRODUCTS OF FERMENTATION OF THE S AND R 
FORMS OF YEASTS! 


By F. W. Fastan, research associate in bacteriology and hygiene, and L. J. WicKEeR- 
HAM, graduate assistant in bacteriology, Section of Bacteriology and Hygiene, 
Michigan Agricultural Experiment Station 


INTRODUCTION 


[In the course of a study with dissociated forms of yeasts it was 
noticed that during the process of fermentation certain forms pro- 
duced a pronounced odor of esters in the fermenting solutions. The 
question naturally arose as to whether the products of fermentation 
of the R forms differed materially either qualitatively or quantitatively 
from the products of fermentation of the corresponding S forms. An 
experiment was therefore set up to determine the amounts of certain 
products present at different times during the process. Yeasts which 
had been dissociated previously ? were used in this experiment. 


CULTURES USED 


The cultures used are described below: 

S form of Saccharomyces aceris-sacchari Fabian and Hall.—This 
yeast produces a glistening, moist growth on agar slants and convex 
colonies on agar plates. The growth is mainly at the bottom in 
liquid media, although a thin film forms on the surface after a few 
days’ growth. The cells average 3.5 to 7 microns in diameter, are 
spherical, and are uniform in appearance. 

R form of Saccharomyces aceris-sacchari.—This yeast produces a 
rugose, white, dry growth on agar slants and flat colonies on agar 
plates. It produces abundant growth on the surface of liquid media 
and a moderate growth at the bottom. The cells average 10 to 16 
microns in length and vary in shape, some being round and others 
elongate. The individual cells do not present a uniform appearance 
as do the S cells, for the R cells contain one or more droplets and a 
large central vacuole. 

S form of Pichia alcoholophila Klécker—This yeast produces a 
glistening, moist growth on agar slants and convex colonies on agar 
plates. It is a bottom yeast. The cells average 6 microns in diame- 
ter and are spherical or slightly elongate. The S cells occasionally 
show small oil droplets. 

R form of Pichia alcoholophila.— This yeast produces a dull, wrinkled 
growth on agar slants and a raised, rugose colony on agar plates. It 
produces an abundant growth both on the surface and at the bottom 
of liquid media. The cells average 10 to 20 microns in length. They 
contain large vacuoles and oil droplets. 





Received for publication Apr. 27, 1936; issued March! 1937. Some of the data in this paper were 
taken from a thesis by the junior author presented in partial fulfillment of the requirements for the degree 
of master oi science, Michigan State College. Journal article no. 244 (N. 8.). 

2 WickeRHAM, L. J., and FABIAN, F. W. DISSOCIATION OF SACCHAROMYCES, FABIAN AND HALL, and 
PICHIA ALCOHOLOPHILA KLOCKER. Jour. Infect. Diseases 58: 165-171. 1936. 
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S form of Saccharomyces cerevisiae Hansen, Saaz strain.— This 
yeast, which for the sake of brevity is referred to in this paper as 
Saccharomyces cerevisiae Saaz, is an industrial strain that produces a 
bottom alcoholic fermentation. The colonies on agar plates are con- 
vex, and agar slant cultures have a moist, glistening appearance 
The cells are spherical and average 6 microns in diameter. 

R form of Saccharomyces cerevisiae Saaz.—This yeast produces a 
wrinkled, white growth on agar slants and a raised, rugose colony. 
The growth is mainly at the surface of liquid media, although there is 
moderate growth at the bottom. The cells are 14 to 20 microns in 
length and vary considerably in shape. Like other R forms, they 
contain oil droplets and large vacuoles. 

SR form of Willia anomala Saito.—The S form of this yeast was 
first obtained by Fabian and McCullough * from the naturally occur- 
ring R form in 1933. It was put away in a stock collection, and when 
it was cultured again for use in this experiment it showed partial 
reversion to the R form. It produces a colony which is Jess wrinkled 
than the true R, is gray instead of white, and appears slightly moist. 
It produces a thin scum on liquid media. Because it has characteris- 
tics of both types it is designated the SR form. 

R form of Willia anomala.—This yeast produces a white, rugose 
growth on agar slants. The colonies are rugose and raised. The 
cells range in length from 14 to 20 microns. The growth in liquid 
media is abundant at the surface and moderate at the bottom. The 
cells show great variation in shape. Only a small percentage of 
them contain oil droplets, but all contain large vacuoles. 


MEDIA USED 


Three media were used: Cider, malt-extract broth, and a synthetic 
medium. The cider had an initial Balling reading of 12°, which was 
increased to 22° when sucrose was added in order to supply sufficient 
sugar for a maximum fermentation by all the yeasts. Dextrose 
instead of sucrose was added to the cider in which the S and R forms 
of Pichia alcoholophila were to be grown, since this yeast ferments 
dextrose very slightly but does not ferment sucrose at all. After 
autoclaving, the cider used for fermentation by the S and R forms of 
P. alcoholophila had a pH value of 3.6, and 0.357 g per 100 cc of total 
acid calculated as malic acid. The cider used for fermentation by the 
two forms of Saccharomyces cerevisiae Saaz had an initial pH value 
of 3.6 and contained 0.400 g of acid per 100 ce calculated as malic. 
The cider used for fermentation by the two forms of S. aceris-sacchari 
and the two forms of Willia anomala had an initial pH value of 3.8 
and contained 0.340 g of acid calculated as malic per 100 cc. 

After the cider had been sterilized but before it was inoculated 
with the yeasts, an analysis was made to determine the presence of 
esters and alcohol. The results were negative. This procedure was 
followed with all the other media used, and with these also the results 
were negative. 

The malt-extract broth consisted of 0.5 percent of sodium chloride, 
0.3 percent of meat extract (Liebig’s), 1.0 percent of peptone (Bacto), 
10 percent of dextrose, and 2 percent of Trommer’s malt extract. 


3 FaBIAN, F. W., and McCuLiouGn, N. B, DISSOCIATION IN YEASTS. Jour. Bact. 27: 583-622, illus. 
1934. 
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Trommer’s malt extract was used because it is high in total solids 
and produces excellent growth. All of the malt-extract broth needed 
was prepared at one time. After autoclaving it contained 0.139 g of 
total acids per 100 ce of medium calculated as lactic acid. The 
pH value was 5.2. 

The synthetic medium contained 0.10 g each of monobasic potassium 
phosphate, monobasic ammonium phosphate, calcium sulphate, and 
magnesium sulphate, 10.0 ¢ of sucrose, and 90 ce of distilled water. 
The reaction of the medium was adjusted with concentrated sulphuric 
acid. This medium also was made up in one batch. After auto- 
claving it had a pH value of 4.5 and contained 0.105 g of total acids 
per 100 ce of medium, calculated as sulphuric acid. 


INOCULATION OF MEDIA 


The media were placed in 6-liter Erlenmeyer flasks, each flask 
receiving 5 liters of the medium. These flasks were fitted with rubber 
stoppers through which a tube of glass and a siphon assembly passed. 
The tube of glass provided for the interchange of gases and was 
plugged with cotton. The siphon provided for the removal of samples 
of media. The latter had two Hoffman clamps separated by a glass 
tube, thus lessening the danger of contaminating the medium. After 
autoclaving, the rubber stoppers were given a heavy coating of 
paraffin. 

Tbe media in the flasks were incoculated with 2 percent of 24-hour 
cultures of the various yeast forms. Extra precautions against con- 
tamination were observed by placing the flasks inside a steamer 
which had been running just previously for several minutes 


CHEMICAL AND YEAST ANALYSIS 


The methods of analysis used, unless otherwise stated, were those 
of the Association of Official Agricultural Chemists.‘ The analyses 
include pH, total acidity, volatile acids, nonvolatile acids, alcohol, 
and ester determinations. At the same time that the chemical 
analysis was made, the cultures were plated on malt-extract agar in 
order to check the cultures for contamination and to see whether any 
dissociation was taking place within the fermenting media, and in 
case there was, to determine the extent of the dissociation. The pH 
was determined by a quinhydrone electrode. Preparation of the 
sample was carried out according to the method of the Association 
of Official Agricultural Chemists, XVIT, 23. The method used for 
determining the total acidity is that given under XVII, 44, except 
that 10 cc of the sample was diluted with 90 cc of distilled water 
before titrating in order to give a clearer end point. Factors of 
0.0067, 0.0049, and 0.0090 were used in calculating the total acidity 
as malic, sulphuric, and lactic, respectively. The volatile acids were 
determined according to the procedure outlined under XVII, 45. 
The nonvolatile acids were determined according to the procedure 
given under XVII, 47, except that they were expressed as malic, 
sulphuric, and lactic acid and factors of 1.116, 0.817, and 1.500, 
respectively, were used. Alcohol was determined by neutralizing 
100 ce of the sample with 20-percent sodium hydroxide solution, dis- 


4 ASSOCIATION OF OFFICIAL AGRICULTURAL CHEMISTS. OFFICIAL AND TENTATIVE METHODS OF ANALYSIS... 
Ed. 3, 593 pp., illus. Washington, D.C. 1930. 
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tilling, and then obtaining the percentage by volume of alcohol by a 


Tag ebullioscope. 


cedure given under XVII, 65. 


from 25 


esters in the sample. 


All analyses were run in duplicate. 


The ester determination is essentially the pro- 
However, it was necessary to use 
25 to 150 ce of excess alkali, depending upon the amount of 


TABLE 1.—Analyses (percent) of cider, malt-extract broth, and a synthetic medium 
containing 10 percent of sucrose or dextrose, during fermentation by the S and R 


forms of Saccharomyces aceris-sacchari, Willia anomala 
J y ’ , 


and S. cerevisiae Saaz 


SYNTHETIC MEDIUM CONTAINING 10 PERCENT SUCROSE, 
8. ACERIS-SACCHARI 


Number of days 


fermented pH 

0 4.5 
15 3.3 
30 2.8 
45 2.6 
60 2.6 
75 2.6 
i) 2.5 

SYNTHETIC MEDIL 
0 4.5 
15 2.7 
30 2.6 
5 2.5 
60 2.5 
75 2.5 
90 2.5 


CIDER CONTAINING 10 





Titrata- | Volatile Non- 
ble acid as | volatile 
acidity acetic acid 
0. 105 0.001 0. 104 
. 176 . 037 . 146 
. 251 . 075 . 190 
. 263 . 128 . 158 
.317 . 198 . 155 
. 343 . 226 . 158 
. 361 . 239 . 166 


1M CONTAINING 10 PERCENT SU 


Alcohol 
by 
volume 


or OO 
x 


8. ACERIS-SACCHARI! 





0. 001 0. 104 
- 029 | . 143 
. 089 . 153 
- 151 . 152 
. 198 . 160 
. 254 . 153 
- 407 . 230 


Newry, 
oe 


Esters 
as ethyl 
acetate 


0. 00 | 


8S FORM OF 


Pichia alccholophila, 


Percentage of S and 
R or SR cells in 


medium 


.03 | 1008. 


05 
-10 

27 
51 
.70 


0. 00 
. 05 
-ll 


27 
43 


. 60 


lw Ss 
100 5 
100 S 
100 8. 
100 8. 


CROSE, R FORM OF 


100 R. 
100 R. 
100 R 


| 100R 


100 R. 
100 R. 





‘ENT SUCROSE, S FORM OF S. ACERIS-SACCHARI? 


0. 3.8 
5 4.0 
15 3.7 
25 3.7 
35 3.6 
45 3.8 
55 3. 6 
65 R.7 
75 BT 


{ 
CIDER CONTAINING 10 


PERC 
0.340} 0.008} 0.331 | 
351 .017| =. 332 
.469| = _.079 | 381 
513 . 129 . 369 
- 521 . 146 358 
. 501 . 158 325 
- 481 145 320 

| .450]) .097 342 

| .981) .0%| .208 


PERCENT SUCROSE, R FORM OF S. ACERIS-SACCHARI? 





0. 008 | 
-017 
042 


0. 331 
. 354 
. 382 
. 409 
. 392 
. 351 
. 344 
. 337 
. 330 


10 PERCENT DEXTROSE, 





> OO 


“Is 
S: 


0. 00 


POW wwe, 
x 
tw 


8. ACERIS-SACCHARI3 


0 3.8 
5 3.8 
15 3.7 
25 3.7 
35 3.6 
45 3.7 
55 3.5 
65 3.5 
75 3.7 
MALT-EXTRACT BROTH CONTAINING 
0 5.2 0 
15 4.1 ; 
30 3.8 ‘ 
45 3.8 ° 
60 peteeadins Geliomantiles | 4.1 
ebateakes aa 3.9 : 


See footnotes at end of table. 





139 0. 002 
383 . 075 
504 . 128 
601 . 238 
. 43 - 224 
371 . 121 


136 
271 
. 312 
. 244 
. 207 
. 190 


0. 


ww 
23S 


on on cn 
+45 


0. 00 
.O1 


87 


0. 00 
00 
Ol 
. 04 
. 06 


0. 00 
. 02 
-10 


- 96 
1.37 


100 8. 
100 5. 
100 8. 
1008 
1S. 
1008 
1008 
1008 


100 R 
100 R 
100 R 
100 R 
100 R 
100 R. 
100 R. 
100 R 


S FORM OF 


958, 91-SR 
78, 93 SR. 
685, 
208, 


63 3, 


42SR 
80 SR. 
73 SR. 























T 


30 
45 


60 


vO 





ore OO DO 


ABLE 1. 


and S. cerevisiae Saaz—Continued 


MALT-EXTRACT 


j 
| Titrata- 


Number of days 
= oH ble 
fermented I acidity 
. 5.2 0. 139 
4.1 . 318 
3.7 578 
3.6 794 | 
3.7 788 
3.7 715 


SYNTHETIC MEDIUM CONTAINING 10 PERCENT SUCROSE, § 


4.5 0. 105 | 
2.8 . 158 
2.9 . 212 
2.5 . 278 
2.4 364 
2.6 . 337 
2.4 - 428 


SYNTHETIC MEDIUM CONTAINING 10 PERCENT SUCROSE, R FORM OF 
ANOMALA! 


4.5 0. 105 
2.8 . 152 
2.9 218 
2.5 271 
2.3 369 
2.4 355 
2.3 443 


CIDER CONTAINING 10 PERCENT SUCROSE, SR! FORM 


3.8 0. 340 
3.9 | . 362 
3. 6 415 
3.7 . 443 
3.6 - 476 
3.6 . 485 
3.5 . 499 
3.5 . 493 
3.7 . 444 


BROTH CONTAINING 
8. ACERIS-SACCHARI3 


Volatile 
acid as 
acetic 


0. 002 
035 
. 088 
- 176 
.214 
- 156 


10 PERCENT 
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Analyses (percent) of cider, malt-extract broth, and a synthetic medium 
containing 10 percent of sucrose or dextrose, during fermentation by the S and R 
forms of Saccharomyces aceris-sacchari, Willia anomala, Pichia alcoholophila, 


DEXTROSE 


Non- Alcohol | Esters 
volatile by as ethyl 
acid volume | acetate 
0. 136 0.00 | 0.00 

. 266 1,12 . 02 

. 347 2.75 | . 04 

. 530 4. 16 | -16 

. 467 4.80 - 42 

. 481 5 .8l 


W. ANOMALA! 


0.001 | 
022 
. 068 
. 169 
291 
281 
. 368 


Ww 


0. 001 
. 025 
. 083 
. 178 
. 288 


0. 008 
010 | 
. 027 


O48 


0. 104 0.00 0. 00 
140 - 50 . 03 

. 157 1.12 .09 

. 140 1. 33 . 25 

. 126 1,19 . 36 

. 107 1. 40 74 

. 127 . 82 . 67 


0. 104 0.00 0.00 
. 132 49 | 03 
. 150 76 | ro 
. 126 . 98 25 
. 134 .76 .29 

. 140 1.12 . 46 | 
namaeai 94 42 
OF W.A 
0. 331 0.00 0.00 
. 351 25 00 
. B85 95 05 
. 389 1. 57 .10 
. 378 2. 08 14 
. 343 | 1. 26 . 28 
2. 22 . 52 
1.80 | . 82 
. 96 





2. 52 | 


CIDER CONTAINING 10 PERCENT SUCROSE, R FORM OF W. 


3.8 0. 340 
4.0 . 363 
3.6 . 423 
a7 . 452 
3.6 510 
3.5 515 
3.5 561 
3. § 545 
3.7 . 444 | 


MALT-EXTRACT BROTH CONTAINING 10 PERCENT DEXTROSE, 
ANOMALA 3 


5. 2 0. 139 

™ 3.5 . 468 
3.5 - 602 

3.5 - 681 

3.7 . 674 

° 3.6 715 





See footnotes at end of table. 


0. 008 
.010 
. 038 
. 075 

107 
. 142 
. 185 
. 184 


. 125 


w. 


0.00 | 0.00 
18 00 
50 . 04 
81 . 09 
1.00 18 
29 25 
1.10 . 40 
47 

1. 22 74 








0. 002 | 
. 040 
. 076 
, ee 
. 174 


0. 136 0.00 | 0. 00 
. 408 1.61 | . 03 
. 488 3. 92 10 
ae 5. 82 .21 
413 6. 35 49 
. 444 5.45 . 88 
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» R FORM OF 


Percentage of S and 
R or SR cells in 
medium 


100 R 
98 R,25R. 
25 R, 75 SR. | 
3 R, 97 SR. 
100 8. 


R‘ FORM OF 


100 SR. 
100 SR. 
100 SR. 
100 SR. 
100 SR. 
100 SR 


| 100 R. 

100 R. 
| 100 R. 
| 100 R. 
100 R. 
100 R. 


NOMALA ? 





100 SR. 

100 SR. 

100 SR. 

100 SR. 

2 R, 80 SR, 18S. 
2 R, 60 SR, 38 S. 
40 SR, 60S. 

7 SR, 93 S. 


ANOMALA ? 


100 R. 
100 R. 
100 R. 
100 R. 
100 R. 
100 R. 
100 R. 
100 R. 


SR‘ FORM OF 


100 SR. 

28,95 SR,3R. 
295, 71 SR. 

73 8, 27 SR. 

60 R, 40 SR. 
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forms of Saccharomyces aceris-sacchari, Willia anomala, Pichia alcoholophila, 


and S. cerevisiae Saaz—Contin 


MALT-EXTRACT BROTH CONT 





1ed 


AINING 


10 PERCENT SUCROSE, 


R FORM OF 








W. ANOMALA 3 
: 4 Titrata- | Volatile Non- Alcohol Esters | Percentaze of 8 and 
- — . . ' pH ble acid as | volatile by asethyl | R or SR cells in 
aimee acidity acetic acid volume acetate medium 
* . 5.2 0. 139 0. 002 0. 136 0. 00 0. 00 
15 4.3 . 249 . 028 . 207 1.04 .05 | 100 R. 
30 3.5 . 567 . 070 . 464 1. 97 .06 | 100 R. 
astade 3.5 . 668 . 163 . 424 2. 50 .23 | 100 R. 
60_. 3.7 - 725 . 186 . 446 3. 02 .96 | 96 R,4SR. 
, Tae hoceindiainnta 3.6 615 . 153 . 386 1.35 1.13 | 99 R,1S8. 
—_— ——— ! : — 
CIDER CONTAINING 10 PERCENT DEXTROSE, S FORM OF P. ALCOHOLOPHILA ? 
0. : 3.6 0. 357 0. 006 0. 351 | 0.00 0.00 | 
5 3.9 . 362 . 008 | . 354 ell | .00 | 1008. 
15 4.1] . 397 . 009 . 338 .09 | .00 | 98S, 2R. 
25 " 3.8 . 344 012 . 331 12 -00 | 768, 24 R. 
35 3.6 . 355 015 . 337 Ki .00 | 268, 74 R. 
45 3.7 . 338 .017 . 319 .18 .00 | 128,88 R 
55 3.9 . 290 014 - 279 | .18 .0 | 10S, 90 R. 
65 3.7 . 239 . 025 - 211 | 12 .00| 78,94R 
75 | 3.9 . 186 .018 . 165 8 .00 | 48,96 R. 


CIDER CONTAINING 10 PERCEN 


ce 
0. 357 


0 3.6 

5 3.8 . 354 
15 3.9 . 351 
25 3.9 . 348 
35 3.6 . 354 
45 3.7 . 337 
55 3.9 . 301 
65 3.7 . 248 
75 3.9 . 217 





0 5.2 
15 5. 1 
30. 1.9 
45 4.8 
60 4.5 
75 4.2 


0 5.2 0. 139 
15 5.1 . 149 
30 4.9 . 207 
45_. 4.8 . 246 
60 4.5 255 
75 4.0 288 


| | 
0 3.6 | 0. 400 
5 3.8 - 480 
15 3.6 . 547 
25 3.6 . 496 
35 3.6 . 515 | 
45 3.6 . 550 
55 3.5 . 493 
65 3.5 | . 484 | 
75 3.3 . 480 | 








351 | 


0. 006 0 0. 00 
. 009 . 344 .10 
. 008 . 342 . 08 
.013 . 333 07 
. 012 . 342 .12 
.014 . 322 .10 
. 016 . 283 .18 
.017 . 229 .31 
. 019 . 195 | .40 


HOLOPHILA 


0. 002 
. 003 
. 009 
.012 
. 008 
. 007 


HOLOPHILA 


0. 002 
. 002 


0 


a 
014 
013 





MALT-EXTRACT BROTH CONTAINING 10 PERCENT DEXTROSE, § 


. 139 0. 00 
. 136 . 20 
. 178 14 
231 - 25 
219 - 22 
- 245 . 20 


. 136 0. 00 
5 | 13 

ehalall . 25 

. 207 21 | 

234 -20 | 

269 25 


T DEXTROSE, R FORM OF P. ALCOHOLOPHILA # 





0. 00 H 
.00 | 100 R. 4 
.00 | 100 R. : 
.00 | 100 R, 

.00 | 100 R, 
.00 | 100 R. 
.00 | 100 R, 

.00 | 100 R, 

.00 | 100 R. 


FORM OF P. ALCO- 


. 00 

-00 | 1008. 

.00 | 95SR,585 

.01 | 72R,148R,148 
.O1 | 100 R. 

.02 | 89SR,7R,48. 


MALT-EXTRACT BROTH CONTAINING 10 PERCENT DEXTROSE, R FORM OF P. ALCO- 






0.00 | 100 R. 
.00 | 100 R. 
.00 | 100 R. 
01 | 100 R. 
Ol | 100 R. 
.00 | 100 R. 











0. 006 0. 0 
. 009 .470 : 
. 046 . 996 1.4 
. 073 415 4. 25 
070 . 437 | 4. 45 
. 080 | . 461 5. 08 
. 077 . 407 | 5. 55 
. 081 395 | 5. 87 
. 069 - 402 | 5.41 


SUCROSE, 


8 FORM OF 8. CEREVISIAE SAAZ? 











393 | 





0. 
































. 00 

. 00 100 S$ 
----| 1008 
.01 100 8 
01 | 1008 
.Ol 1008 
.O1 (5) 
.00 | (6) 
-00 | (®). 











See footnotes at end of table 
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-Analyses (percent) of cider, malt-extract broth, and a synthetic medium 


containing 10 percent of sucrose or dextrose, during fermentation by the S and R 
forms of Saccharomyces, aceris-sacchari, Willia anomala, Picha alcoholophila, 
and 8. cerevisiae Saaz—Continued. 


CIDER CONTAINING 10 PERCENT SUCROSE, R FORM OF 8S. CEREVISIAE SAAZ? 


- “a Titrata- | Volatile Non- Alcohol Esters | Percentage of Sand 
oe oe pH ble acid as | volatile by as ethyl | R or SR cells in 
. acidity acetic acid volume | acetate medium 
3.6 0. 400 0. 006 0. 393 0. 00 0. 00 
3.9 473 012 | . 460 31 .00 | 100 R. 
3.5 . 029 | . 456 -77 .03 | 100 R. 
25 3.6 . 047 - 440 1, 25 -07 | 100 R. 
35 3.6 . 070 . 446 | 1. 28 .14 | 100 R. 
5 3.5 .112 | . 454 | 1. 46 .22 | 100 R. 
55 3.3 - 128 | . 448 1.90 .36 | 100 R. 
65 3.3 . 105 | . 404 1, 21 -52 | 100 R. 
75 3.1 . 064 . 408 . 78 1.00 | 100 R. 











MALT-EXTRACT BROTH CONTAINING 10 PERCENT DEXTROSE, S FORM OF § 


CEREVISIAE SAAZ 


3 


0 5.2 0. 139 0. 002 0. 136 | 0.00 0. 00 

15 4.3 - 292 . 047 . 222 4. 66 .00 | 1008. 

40) 4.3 . 341 . 054 . 260 6. 25 .00 | 1008, 

45 4.3 309 . 054 . 228 6.72 .01 | () 

60 4.3 . 295 . 060 . 205 | 6. 26 -01 | 95S, 5 R.? 
75 4.1 . 296 - 059 | . 207 | 6.72 .02 | 74SR, 268.8 





MALT-EXTRACT BROTH CONTAINING 10 PERCENT DEXTROSE, R FORM OF § 


CEREVISIAE SAAZ3 





) 5.2 0. 139 0. 002 0. 136 0. 00 0.90 | 100 R. 
15 4.4 . 249 . 030 204 1. 48 .02 | 100 R. 
0) 4.1 . 374 . 072 . 266 3.00 09 | 100 R. 
45 1.0 . 476 067 .217 4. 35 .27 | 100 R. 
er 3.8 . 585 . 193 . 296 4. 51 .72 | 100 R. 
75 = 3.9| 2504 "162 261 5.18 1:07 | 100 R. 








Titratable acidity and nonvolatile acid as sulphuric. 
2 Titratable acidity and nonvolatile acid as malic. 
’ Titratable acidity and nonvolatile acid as lactic. 
‘Some revision of S to R Form. 


5 No growth in 1 loopful. 

* No growth in 5 ce 

’ Seant growth in 5 ce. 

8 Moderate growth in 5 ec, 


RESULTS OF EXPERIMENTS 
ANALYSES 


The data given in table 1 show that in the cider and in the malt- 
extract broth the S form of Saccharomyces aceris-sacchari reached a 
peak in volatile acid production by the forty-fifth day. This peak 
was followed by a rapid decline in the concentration of the volatile 
acids until the seventy-fifth day, when approximately 50 percent of 
the maximum amount remained. In the synthetic medium (table 1) 
the volatile acids increased steadily up to the ninetieth day, however, 
the concentration at that time equaling the concentration of the vol- 
atile acids in the malt-extract broth when at their peak. It is inter- 
esting to note that, although the maximum production of volatile 
acids in the cider reached little more than half the maximum produc- 
tion in either of the other two media, a decline in concentration of the 
volatile acids nevertheless took place. Table 1 shows that by the 
seventy-fifth day the percentage of alcohol was practically the same 
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in the synthetic medium and in malt-extract broth, but that the cider 
contained between 40 and 50 percent more alcohol. Ester production 
was 24 percent greater in the cider than in the synthetic medium, and 
96 percent greater in the malt-extract broth than in the synthetic 
medium. 

During the 75-day fermentation period dissociation occurred in the 
yeasts. However, the increase in ester production is directly attrib- 
utable to the action of the S forms on the malt-extract medium rather 
than to the intermediate forms of yeasts, since, as will be shown later, 
the R form of this yeast produced less esters than did the S form. 

Figure 1, A, shows the relationship between alcohol, volatile acids, 
and esters produced by the S form of Saccharomyces aceris-sacchari in 
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FIGURE 1.—Relationship between the alcohol, volatile acids, and esters produced by the S and R forms 
of Saccharomyces aceris-sacchari and S. cerevisiae Saaz in cider; the ratios used in plotting alcohol, esters, 
and acids being 1:2:20, respectively; A, S form of S. aceris-sacchari; B, R form of S. aceris-sacchari; C, § 
form of S. cerevisiae Saaz; D, R form of S. cerevisiae Saaz. 


cider. It will be noted that the amount of volatile acids and alcohol 
increased rapidly after the fifth day, though the ester production was 
very scant until approximately the fortieth day, after which it 
increased rapidly. 

Table 1 shows that the volatile acids produced by the R form of 
Saccharomyces aceris-sacchari, as well as by the S form, increased con- 
tinuously in the synthetic medium throughout the experiment, but 
decreased after 50 or more days in cider and malt-extract broth. 
Unlike the condition obtaining with the S form, the concentration of 
volatile acids produced by the R form in the synthetic medium at the 
seventy-fifth day was about 6 percent greater than the greatest con- 
centration of this product in either cider or malt-extract broth. At 
60 days the alcohol production had reached its peak in all three media. 














jan. 15,1937 Products of Fermentation oi Sand R Forms of Yeasts 155 











The cider (shown at 65 pains showed a ee alcohol concentration of 
5.60 percent, malt-extract broth 4.80 percent, and synthetic medium, 
3.47 percent by volume. The ester production in all three media was 
about the same. Figure 1, B, shows graphically the production of 
fermentation products by the R form in cider. It will be noted that 
about 15 days after the production of esters started there was a rapid 
decrease in the volatile acids and later a decrease in the alcohol. 

The amount of volatile acids produced by the SR form of Willia 
anomala (table 1), was approximately equal in the cider and in the 
malt-extract broth at the seventy-fifth day, but in the synthetic 
medium the amount produced was 61 and 65 percent more than in 
the cider and malt-extract broth, respectively. Alcohol production 
varied greatly in the different media. At the seventy-fifth day the 
amount of alcohol present in the synthetic medium was 1.40 percent, 
in cider 2.52 percent, and in the malt-extract broth 5.45 percent. The 
relatively large amount of alcohol in the last-mentioned medium may 
have been due in part to the presence of a large percentage of true S 
cells which appeared in the medium during the latter part of the fer- 
mentation. The difference in ester production was not large. The 
cider had the highest concentration of esters and the synthetic medium 
the lowest. 

The R form of Willia anomala reached a peak of 0.174 and 0.186 ¢ 
of volatile acids per 100 ce in malt-extract broth and cider, respec- 
tively, at the sixtieth day (table 1). In synthetic medium a peak of 
().288 g per 100 cc was reached at the same time. In the malt-extract 
broth ‘approximately 135 percent more alcohol was produced than in 
the cider and 170 percent more than in the synthetic medium. These 
differences can only be due to the composition of the media, as prac- 
tically no dissociation took place in any of the three media. The 
amount of ester production on the seventy-fifth day was 0.46, 0.74, 
and 1.13 g per 100 cc in the synthetic medium, cider, and malt-extract 
broth, respectively. 

None of the S and R forms of Pichia alcoholophila or Saccharomyces 
cerevisiae Saaz grew in the synthetic medium. Therefore, only the 
data concerning the fermentation of these yeasts in cider and malt- 
extract broth are presented. 

Table 1 indicates the slight fermentation produced by the S form 
of Pichia alcoholophila. In malt-extract broth the greatest concen- 
tration of volatile acids was 0.012 g per 100 cc, while in cider it was 
0.025 g per 100 cc. The alcohol concentration in both media varied 
considerably during the fermentation, but did not rise above 0.48 
percent. Only traces of esters were produced. 

The fermentation produced by the R form of Pichia alcoholophila 
(table 1) was nearly the same as the fermentation produced by the S 
form, but more volatile acids were produced in cider than in the malt- 
extract broth. The slight fermentations produced by the S and R 
forms were not due to poor growth of the yeast forms, since both 
grew abundantly in the two media used. 

Table 1 shows that the maximum amount of volatile acids produced 
by the S form of Saccharomyces cerevisiae Saaz in the cider exceeded 
that produced in malt-extract broth by 35 percent. The alcohol pro- 
duction was much more rapid in malt-extract broth and reached a 
higher maximum than in cider. Only a very small amount of esters 
was produced, the highest concentration being 0.02 g per 100 cc in 
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malt-extract broth. When one loopful of malt-extract broth was 
plated on the forty-fifth day, no colonies resulted. At 60 days 5 cc of 
this medium produced only a few colonies when plated. At the 
seventy-fifth day 5 ce of the culture produced about 400 colonies. 
The culture was dead in cider by the fifty-fifth day. Plates inoculated 
with 5 ce of the medium produced no colonies. Figure 1, C, shows 
the course of fermentation by the S form in cider. 

Table 1 shows that the R form of Saccharomyces cererisiae Saaz 
produced 50 percent more volatile acids and 173 percent more alcohol 
in malt-extract broth than in cider. The ester production was nearly 
the same in the two media, about 1 g per 100 cc. Figure 1, D, shows a 
rapid decrease in alcohols and volatile acids after the fifty-seventh 
day. The esters, however, continued the rapid rise which started 
about the twenty-fifth day. 

The last column in table 1 shows the amount of dissociation which 
took place in terms of the percentage of cells that had reverted at the 
end of 60 days. In no case was there any change of form in the 
synthetic medium. Four of the forms which dissociated in malt- 
extract broth did not dissociate in cider. Furthermore, the malt- 
extract broth contained a greater number of dissociated forms at 60 
days than did the cider in which the same species and form of yeast 
Was growing. 

INFLUENCE OF OXYGEN ON ESTER PRODUCTION 


In order to study the influence of oxygen on ester production an 
experiment was set up in which the amount of air in contact with the 
fermenting solution varied considerably. For this experiment the R 
form of Saccharomyces cerevisiae Saaz was seeded in malt-extract 
broth containing 10 percent of dextrose and 2 percent of malt extract. 
Previous experiments had demonstrated that this combination of 
yeast and medium would produce a high ester content. A 6-liter 
flask containing 5 liters of the medium was seeded with the yeast. 
The only access of air was through a hole 15 mm in diameter which 
was tightly stoppered with cotton. The flask was allowed to remain 
undisturbed during the fermentation. At the end of 10 days less than 
0.02 g per 100 ce of esters had formed. Five liters of medium was 
placed in a second 6-liter flask and the 50-mm opening to the flask 
plugged with cotton. After inoculation and every other day thereafter 
the flask was shaken in order to incorporate air into the medium. 
At the end of 10 days 0.58 g of esters per 100 cc had been produced. 

A third 5 liters of medium was placed in a sterile pan 16 inches 
square. The medium covered the pan to a depth of about 1 inch. 
To prevent contamination the pan was covered with several layers of 
sterile cheesecloth. The contents were not agitated during the fer- 
mentation. The yeast made an extremely abundant growth, but at 
the end of 10 days there were no esters present. 

This experiment indicates that the conditions under which yeasts 
are grown influence the amount of esters produced. When oxygen is 
low, as in the first flask, ester production is low. When oxygen is 
abundant, as in the pan, no esters are formed. The ideal condition, 
or at least the most favorable condition of the three tried here, was 
when the amount of oxygen present was sufficient to maintain a 
normal growth, as in the second flask. 
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DISCUSSION 







The comparison of the products of fermentation of the S and R 
forms of Saccharomyces aceris-sacchari as given in table 1, shows that 
the medium is a more important factor than the form of yeast in 
determining the amount of volatile acids produced. To illustrate, 
the R form produced much more volatile acids in the synthetic 
medium and in cider than the S form did, but the S form produced 
slightly more in the malt-extract broth than did the R form. Without 
exception the S form produced more alcohol than did the R form in 
the same kinds of media. In cider and malt-extract media the S 
form produced more esters than the R form; the reverse was true in 
the synthetic medium. 

Table 1 shows no significant difference in the production of vola- 
tile acids in the three media between the SR and R forms of Willia 
anomala. In all media the SR form produced more alcohol than 
the R form. The difference was not great in the synthetic medium, 
but in malt-extract broth the SR form produced 110 percent more 
alcohol than did the R form. Ester production varied with the kind 
of medium used. In malt-extract broth the R form produced 28 
percent more esters than did the SR form, but the SR form pro- 
duced 61 percent more in the synthetic medium and 30 percent more 
in the cider. 

Table 1 indicates that there was no significant difference between 
the products of fermentation of the S and R forms of Pichia 
alcoholophila. 

The R form of Saccharomyces cerevisiae Saaz (table 1) produced 
221 percent more volatile acids in malt-extract broth and 58 per- 
cent more in cider than did the S form. The S form far exceeded 
the R form in alcohol production in both media, the ratio in cider 
being approximately 3:1 and in malt-extract broth 1.3:1.0. The S 
form produced practically no esters, but the R form produced ap- 
proximately 1 ¢ per 100 cc in both media. 

The S and SR forms of the four yeasts vary to a marked degree 
in their fermentations. Table 1 shows that the ratio of the maxi- 
mum production of volatile acids in malt-extract broth by the S 
or SR forms of Saccharomyces aceris-sacchari, Willia anomala, S. 
cerevisiae Saaz, and Pichia alcoholophila is 20:15:5:1, respectively. 

The ratio for maximum alcohol production in the same order is 
23:25:27:1, and the maximum ester ratio is 68:44:1:1. Figure 1, 
B, shows the fermentation of cider produced by the R form of Sac- 
charomyces aceris-sacchari, which is a strong ester producer. It will 
be noted that about 15 days after the formation of esters began, the 
volatile acid concentration dropped rapidly. Figure 1, C, shows the 
fermentation of the same medium by the S form of Saccharomyces 
cerevisiae Saaz, which is a very weak ester producer. The volatile 
acids after 25 days of fermentation reached a rather constant con- 
centration which was maintained throughout the rest of the fer- 
mentation. The drop in concentration of the volatile acids in the 
case of the S form of S. aceris-sacchari may have been due to the 
combination of volatile acids and alcohol to form esters. The drop 
in the alcohol concentration would not be as noticeable as the drop 
in volatile acids concentration because of the comparatively large 
amount of the former in the media. B and D, of figure 1, which 
are graphs of fermentations caused by ester-producing yeast forms, 
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show the drop in volatile acid and alcohol which followed the rapid 
increase in ester production. 

Table 1 shows that the R forms differ as much among themselves 
as do the S forms. The ratio of the maximum production of volatile 
acids in malt-extract broth by the R forms of Saccharomyces aceris- 
sacchari, Willia anomala, S. cerevisiae Saaz, and Pichia alcoholophila, 
in the order named, is 8:7:7:1; of alcohol, 19:12:21:1; and of esters, 
81:113:107:1. 

SUMMARY 


Detailed analyses of the products of fermentation of the S and 
R forms of the four yeasts studied, Saccharomyces cerevisiae Saaz, 
S. aceris-sacchari, Pichia alcoholophila, and Willia anomala in three 
different media, cider, malt-extract broth, and a synthetic medium 
in general show: 

(1) In some species of yeasts there is a marked difference in the 
fermentation properties of the S and R forms; in other species there 
is little or no difference. 

(2) The S forms produced more alcohol in a shorter length of 
time than did the R forms. 

(3) The production of volatile acids was extremely variable. 
The amount produced varied with the different media, with the 
species of yeast, and with the S and R forms of the different species. 

(4) The production of esters by the S and R forms of the four 
species of yeasts does not show as great variation as the productions 
of volatile acid. There are some exceptions to this. The most 
important is the S form of S. cerevisiae Saaz in malt-extract broth 
and cider. In these media the amount of esters produced by the 
S form is practically nil while the amount produced by the R form 
is 1 g per 100 ec. Another exception is in malt-extract broth where 
the S form of S. aceris-sacchari produced 69 percent more esters 
than the R form. The SR form of W. anomala produced 30 percent 
more esters in cider than the R form and 61 percent more in the 
synthetic medium. P. alcoholophila produced practically no esters 
in any of the media. 

(5) The conditions under which the yeasts are grown influence the 
ester production. Sufficient oxygen to maintain normal growth pro- 
duced the largest amount of esters. When there was an abundance 
of oxygen or insufficient oxygen no esters were produced. 

(6) In the yeasts which produced esters, the esters did not appear 
in appreciable quantities until 35 to 45 days after inoculation. Shortly 
after the esters appeared there was a corresponding decrease in 
alcohol and volatile acids. Jn the case of the S form of S. cerevisiae 
Saaz there was a moderate volatile acid and high alcoholic content 
in malt-extract broth, but practically no esters were formed. 

(7) The three media used affected the stability of the yeasts dif- 
ferently. The S and R forms of the species which grew in the syn- 
thetic medium were very stable in that medium, the S form showed 
no tendency to change to the R form, nor the R form to change to 
the S form. In the cider, however, the tendency to change was 
greater and in the malt-extract broth it was still greater. 
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